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Jjif ABSTRACT 

O , Full range Herschel/PACS spectroscopy of the (ultra)luminous infrared galaxies NGC 4418 and Arp 220, observed as part of the 

+3 ■ SHINING key programme, reveals high excitation in H2O, OH, HCN, and NH3. In NGC 4418, absorption lines were detected with 

^ ■ £io„ei >800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp 220 the excitation is somewhat lower While outflow 

^ ■ signatures in moderate excitation lines are seen in Arp 220 as have been seen in previous studies, in NGC 4418 the lines tracing its 

outer regions are redshifted relative to the nucleus, suggesting an inflow with M<12 Mq yr"'. Both galaxies have compact and warm 
■ (rdust^ 100 K) nuclear continuum components, together with a more extended and colder component that is much more prominent and 

^ ■ massive in Arp 220. A chemical dichotomy is found in both sources: on the one hand, the nuclear regions have high H2O abundances, 

QQ ■ ~10"^, and high HCN/HtO and HCN/NH3 column density ratios of 0.1-0.4 and 2-5, respectively, indicating a chemistry typical of 

evolved hot cores where grain mantle evaporation has occurred. On the other hand, the high OH abundance, with OH/H2O ratios of 
~0.5, indicates the effects of X-rays and/or cosmic rays. The nuclear media have high surface brightnesses (>10'^^ Lg/kpc^) and are 
estimated to be very thick (A^h^IO""^ cm"-). While NGC 4418 shows weak absorption in H^**0 and '^OH, with a "0-to-'**0 ratio 
of >250-500, the relatively strong absorption of the rare isotopologues in Arp 220 indicates "*0 enhancement, with "'0-to-'*'0 of 
70-130. Further away from the nuclear regions, the H2O abundance decreases to <10 ^ and the OH/H2O ratio is reversed relative 
>— ' ■ to the nuclear region to 2.5-10. Despite the different scales and morphologies of NGC 4418, Arp 220, and Mrk 231, preliminary 

evidence is found for an evolutionary sequence from infall, hot-core like chemistry, and solar oxygen isotope ratio to high velocity 
outflow, disruption of the hot core chemistry and cumulative high mass stellar processing of '*0. 

Key words. Line: formation - Galaxies: ISM - ISM: jets and outflows - Infrared: galaxies - Submillimeter: galaxies 

I 

5" ■ 1 ■ Introduction atoms and ions decrease in strength, as absorption in molecular 

lines, specifically H2O and OH, increase in depth. 

From the first CO observations of IRAS galaxies (Young et 

al., 1984), much eflFort has been devoted to understanding spec- Water (H2O) and hydroxyl (OH) are indeed key molecu- 
troscopic observations of molecules in (ultra)luminous infrared lar species that trace relevant physical and chemical properties 
galaxies (hereafter, (U)LlRGs), and how they are related to the of the interstellar medium in luminous infrared extragalactic 
conditions traced by the continuum and line observations of sources. Powerful OH mega-masers are common in ULlRGs, 
atoms and ions. The launch of the Infrared Space Observatory and H2O masers are found in AGN accretion disks and in post- 
(ISO) and, more recendy, the Herschel Space Observatory, have shocked gas associated with nuclear jets (see Lo, 2005, for a re- 
opened a new window to study extragalactic sources in the view). With many transitions lying at far-lR wavelengths, where 
far-infrared (hereafter, far-lR) domain, where the bulk of the the bulk of the luminosity is emitted, H2O and OH are mainly 
(U)LIRG luminosity is emitted. The relevance of far-lR molec- excited through absorption of far-lR photons (Gonzalez- Alfonso 
ular spectroscopy of (U)LIRGs is well illustrated with the initial et al., 2004, hereafter G-A04), and thus their excitation helps 
ISO result of Fischer et al. ( 1 999), that the fine-structure lines of probe the general properties of the underlying far-lR continuum 

sources (Gonzalez- Alfonso et al., 2008). Chemically, H2O is ex- 

pected to trace an undepleted chemistry where grain mantles are 

* Herschel is an ES A space observatory with science instruments pro- evaporated; both H2O and OH are also boosted in shocks and, 

vided by European-led Principal Investigator consortia and with impor- especially OH, in X-ray and Cosmic Ray Dominated Regions 

tant participation from NASA. (XDRs and CRDRs). Furthermore, OH has turned out to be 
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a unique tracer of massive molecular outflows (Fischer et al., 
2010; Sturm et al., 2011). Recently, H2O has been detected in 
high-z sources (Impellizzeri et al., 2008; Omont et al., 2011; Lis 
et al., 2011; van der Werf et al., 2011; Bradford et al., 2011), 
and thus understanding its role in the local Universe will be cru- 
cial for future routine observations of H2O with ALMA in the 
far Universe. HCN is another key molecule widely studied in 
local (U)LIRGs at (sub)millimeter wavelengths (e.g., Aalto et 
al., 1995). NH3, first detected in extragalactic sources by Martin 
& Ho (1979), is another important tracer; in Arp 220, it was 
detected at far-IR wavelengths with ISO (G-A04) and at cen- 
timeter wavelengths by Takano et al. (2005) and recently by Ott 
et al. (2011). With its high sensitivity and spectral resolution at 
far-IR wavelengths, the Herschel/PACS instrument (Pilbratt et 
al., 2010; Poglitsch et al., 2010) provides new opportunities for 
using the far-IR transitions of these molecules to probe the kine- 
matic, chemical, and radiative activity in IR luminous galaxies. 

We report in this work Herschel/PACS spectroscopic ob- 
servations of H2O, OH, their '^O isotopologues, NH3, and the 
surprising detection of highly excited HCN, in NGC 4418 and 
Arp 220. NGC 4418 is a peculiar, single nucleus galaxy, with a 
moderate luminosity (~ 10" Lq) but with other properties sim- 
ilar to warm ULIRGs: a high Lfir/A^Hj ratio and an extreme [C 
11] deficit (e.g. Gracia-Carpio et al., 201 1), a high continuum sur- 
face brightness, an extremely compact luminosity source (Evans 
et al., 2003), and warm infrared colours {S 2511m/ S eu/jm - 0.23, 
S eotim/ S looi^m - 1.37). We adopt a distance to the source of 29 
Mpc. Arp 220 is the prototypical and nearest ULIRG, showing a 
double nucleus and a luminosity of a; 1.3 x 10'^ L©. We adopt a 
distance to the source of 72 Mpc. The observations are described 
in §2; models are shown in §3; model results are discussed in §4, 
and our main conclusions are summarized in §5. 



2. Observations and results 

2.1. Observations 

The fuU (52.3-98, 104.6-196 yum), high resolution, PACS spec- 
tra of NGC 4418 and Arp 220, taken as part of the guaranteed- 
time key programme SHINING, were observed on July 27th and 
February 27th (2010), respectively. The spectra were taken in 
high spectral sampling density mode using first and second or- 
ders of the grating. The velocity resolution of PACS in first order 
ranges from x 320 km s"' at 105 fim to x 180 kms"' at 190 /jm, 
and in second order from =s 210 km s"' at 52 fim to =s 110 km s ' 
at 98 microns. For NGC 4418, we also present a scan around the 
[O i] 63 i-im line taken in third order with velocity resolution of 
X 85 km s ' . The data reduction was mostly done using the stan- 
dard PACS reduction and calibration pipeline (ipipe) included in 
HIPE 5.0 975. The reduction of Arp 220 required, however, two 
additional steps in order to correct for pointing errors that re- 
sulted in significant variations of the continuum level of the cen- 
tral spaxel' as well as in wavelength off'sets generated by point- 
ing errors along the dispersion direction of the slit (Poglitsch et 
al., 2010). These effects were produced by the movement of the 
source between the central spaxel and a neighbouring spaxel. In 
order to correct for them, (/) the emission from the two spax- 
els with highest signal were added together in an intermediate 
step in the pipeline, which significantly improved the continuum 
baseline of the combined spectrum but increased the noise by a 



' PACS has 5x5 spatial pixels (spaxels), each with 9" x 9" resolution 
and covering a field of 47" x 47" (Poglitsch et al., 2010). 




10^ 10^ 10^ 



Fig. 1. Spectral energy distribution of NGC 4418 and Arp 220 from 
mid-IR to millimeter wavelengths. Both the Herschel/PACS and the 
Spitzer/IRS spectra are shown. Data points at (sub)millimeter wave- 
lengths for NGC 4418 are from (Yang & Phillips, 2007, 350 yum), 
(Roche & Chandler, 1993, 450, 800, and 1100 fim), (Dunne et al., 
2000; Lisenfeld et al., 2000, 850 yum), and (Albrecht et al., 2007, 1300 
fim). For Arp 220, they are taken from (Bales et al., 1989, 450 //m); 
(Rigopoulou et al., 1996, 350, 800, and 1100 yum); and (Sakamoto et 
al., 1999, 2008, 1300 fim). Models discussed in §3 are included. For 
NGC 4418, the yellow, blue, light-blue, and solid-green curves show 
the models for the hot, core, warm, and extended components, respec- 
tively, with parameters given in Table 1 . An additional cold component, 
shown with the dashed-green curve, is included to better fit the SED at 
long wavelengths. For Arp 220, the yellow, blue, light-blue, and solid- 
green curves show the models for the hot, western nucleus, eastern nu- 
cleus, and extended components, respectively (Table 1). 

factor^ a; V2; (//) the time evolution of the pointing shifts was 
reconstructed from the star tracker frames, and wavelength cal- 
ibration was corrected according to Fig. 8 by Poglitsch et al. 
(2010). Pointing errors across the slit were found from to ~ 3", 
resulting in velocity shifts of < 60 kms"'. We estimate that the 
associated uncertainties are below 30 kms"', though they may 
be higher for some individual lines. 

For the final calibration, the spectrum was normalized to the 
telescope flux and recalibrated with a reference telescope spec- 
trum obtained from observations of Neptune. The flux density 
of the spectrum in the central spaxel(s) was then rescaled to the 

- Since the neighbouring spaxel has a signal much lower than the 
central one, the subsequent final calibration leaves the factor V2 nearly 
unchanged. 
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Fig. 2. Energy level diagrams of OH and H2O (ortho and para). Red arrows indicate lines detected in both sources; blue arrows denote lines 
detected only in NGC 4418, green-dashed arrows mark transitions contaminated (blended) by lines of other species, and blue-dashed arrows 
denote marginal detections in NGC 4418. The small numbers in the OH diagram indicate rounded wavelengths in yjm. 



value of the total continuum of the source obtained by the combi- 
nation of the 25 spaxels of the PACS spectrometer. This method 
assumes that the distribution of the IR continuum and the lines 
are similar for the spatial resolution of the observations. Indeed 
the line to continuum ratio in the central and the combined 25 
spaxels was found to be the same for all strong lines except for 
the ground-state OH 119 and 79 fim in Arp 220, which display 
9% and 3% deeper absorption in the central spaxels. The lat- 
ter OH lines thus appear to show low-level reemission (proba- 
bly through scattering) over kpc spatial scales. In general, the 
comparison indicates that both sources are at far-IR wavelengths 
point-sources compared with the PACS 9" x 9" spatial resolu- 
tion (~ 1.3 X 1.3 kpc^ for NGC 4418 and ~ 3.1 x 3.1 kpc^ for 
Arp 220). 

The continua of both sources are shown in Fig 1, including 
the Spitzer IRS spectra (Spoon et al., 2007); the models are dis- 
cussed in §3. Both galaxies show strong silicate absorption at 9.7 
fim, with PAH emission that is weak in Arp 220 and only weakly 
detected at 1 1.2 /urn in NGC 4418 (Spoon et al., 2007). Besides 
the obvious differences in luminosity, a striking contrast between 
both sources is the S (,oiim/S25iim ratio, which is a; 4 in NGC 4418 
and approaches a; 13 in Arp 220. The higher value in Arp 220 
indicates that larger amounts of dust surround the nuclei in an 
extended region (ER in G-A04, denoted here as Cextended), both 
extinguishing the nuclear emission at 25 fim and contributing to 
the far-IR continuum. This emission can be interpreted as heat- 
ing of dust by the radiation emanating from the nuclei (Soifer 
et al., 1999, G-A04). In NGC 4418, the warmer spectral energy 
distribution (SED) suggests that the amount of dust in Cextended is 
lower, but still probably significant (see §3). On the other hand, 
the (sub)millimeter emission in Arp 220 is dominated by the nu- 
clei as revealed by interferometric (sub)millimeter observations 
(Downes & Eckart, 2007; Sakamoto et al., 2008), indicating ex- 



treme opacities in, and high luminosities from the nuclei. No 
interferometric observations are available for NGC 4418. 

Spectroscopic data used for both line identification and ra- 
diative transfer modeling were taken from the JPL (Pickett et 
al., 1998) and CDMS (Muller et al., 2001, 2005) catalogues. 
Figure 2 shows the energy level diagrams of OH and H2O. The 
spectra around the H2O lines are shown in Fig. 3 (blue and green 
bands, A < 100 fj.m) and Fig. 4 (red band, A > 100 fim), while 
high-lying H2O lines tentatively detected in NGC 4418 (blue- 
dashed arrows in Fig 2) are shown in Fig. 5. The OH line pro- 
files are presented in Fig. 6. Figs. 7 and 8 show excerpts of the 
PACS spectra around the Hj^O and "^OH lines, respectively. The 
HCN and NH3 lines are displayed in Figs. 9 and 10, respectively. 
Figure 1 1 compares the line profiles of the [O i] 63 and 145 jum 
transitions with those of some OH lines. Figure 12 compares the 
line shapes of two HCO^ transitions observed with high spatial 
resolution by Sakamoto et al. (2009) in Arp 220, with selected 
lines observed with Herschel/PACS. Results of the continuum 
and Une modeling are summarized in Tables 1 and 2, and the 
line equivalent widths of the H2O, [O i] and OH, HCN, NH3, 
H^**0, and '**OH lines are listed in Tables 3 to 8. 

2.2. Systemic redshifts 

Line parameters were derived by subtracting baselines and fit- 
ting Gaussian curves to all spectral features. One relevant ques- 
tion in our study is the zero-velocity adopted for the high- 
excitation molecular gas in the nuclear region of the galaxies. 
In NGC 4418, the excited H2O lines do not show any evidence 
of systematic shifts indicative of outflowing or inflowing gas. 
Fig 13 indicates that the only H2O line clearly redshifted relative 
to all others is the ground-state 2i2— » loi at 179.5 fim, which is 
expected to be contaminated by CH^ (2-1) (see below). From all 
other H2O lines, the average redshift is z = 0.007035 + 0.000070 



3 



Gonzalez-Alfonso et al.: Herschel/PACS spectroscopy of NGC 4418 and Arp 220 



o 




0.7 



53 



53.1 53.2 



56.2 56.4 57.3 



57.6 



57.9 



0.9 - 



0.7 





-3 ' 






L \ 

- 





4 ' 


-4o4 . - 






r , , 





58.6 58.8 61.5 

1 I 1 r- 



61.8 




1 

:h 331- 


"^20 33Q- 


3o3 








1 









8l8" 


'''07 ^OE 


-7 ' ■ 

17, 










] 








63.2 63.4 63.6 



o 

^ 0.9 



66.6 67 67.2 67.4 



0.1 
0.7 

1 

0.9 
0.8 



;j7i,- 


' 1 ' 1 - 

n ^ n n n n „ ' 






n 




,1 


1 ^ 



71.6 71.8 72 

\ ' I ' I 

m 6ie-5o5 





71.4 



75.6 78.6 78.8 79 



- 


716- 


-"^07 n 















83.4 84.6 84. { 



89.8 



90 



94.6 94.8 



- ' 1 


1 ' - 
-4o4 440-431 : 








~ , 1 


1 




95,6 95.8 



Fig. 3. H2O lines in the blue and green bands {X < 100 lum) detected in NGC 4418 (upper spectra in each panel) and Arp 220 (lower spectra). The 
spectra, histograms with black solid lines, have been scaled by a factor of 0.5 in panels k and p. Vertical dotted lines indicate the rest wavelengths 
of the lines relative to z = 0.00705 and 0.0181 for NGC 4418 and Arp 220, respectively. Model results are also shown. For NGC 4418, the blue and 
light-blue curves show the models for the Ccorc and C„.^ components, respectively; for Arp 220, the blue and light-blue lines show the estimated 
contributions from C„est and Ceast, and the green curve (in panel k) shows the contribution by Ccxtcndcd- Red is total. Model parameters are given in 
Tables 1 and 2. 
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Fig. 4. H2O lines in the red band (A > 100 yum) detected in NGC 4418 (upper spectra in each panel) and Arp 220 (lower spectra). The spectra have 
been scaled by a factor of 2 in panels a,k and 1, and by 3 in panel c. Model results are also shown. Color codes for the contributions of the different 
components are given in Fig. 3, with the yellow curve in panel q showing the additional contribution by the halo component (Chaio) in front of the 
nuclear region of Arp 220. 
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(cz = (2110 + 20) kms"'), where the error denotes the stan- 
dard deviation of the mean. The inferred redshift is fully con- 
sistent with observations of CO, CN, HNC, and other molec- 
ular species observed at millimeter wavelengths (e.g. Aalto et 
al., 2007a). The Sloan Digital Sky Survey (SDSS) specti'um of 
NGC 4418 also gives a consistent redshift for the optical ab- 
sorption/emission lines of z = 0.0070 and 0.0072, respectively 
(D. Rupke, private communication). Finally, this is also consis- 
tent with the redshift measured for the high-lying OH lines (see 
below), which also trace the nuclear component, and with that 
of the [C ii] emission Hne at 157.741 yum, z = 0.007050, which 
we adopt as the reference value. In Arp 220, Sakamoto et al. 
(2009) detected P-Cygni profiles in excited submillimeter lines 
of HCO^ toward both nuclei, indicating dense outflowing gas 
in the nuclei (see also Fig. 12 and §2.7). The transition from 
absorption to emission occurs at a velocity of » 5420 kms"' 
toward both nuclei, which we adopt as the systemic redshift 
(z = 0.0181). The foreground absorption traced by CO (3-2) to- 
ward the western nucleus lies at a lower velocity, 5350 km s ' . 
The line velocity shifts (Vshitt) relative to the above systemic red- 
shifts, the line widths, and the equivalent widths (Weq) are plot- 
ted in Fig. 13 for the H2O lines and in Fig. 14 for the OH lines. 

2.3. H20anclOH 

Based on the H2O spectra and the molecular line positions 
shown in Figs. 3 and 4, most H2O lines are free of contamina- 
tion from lines of other species, with the following exceptions: as 
mentioned above, the H2O 2i2— > loi hne in Fig. 4q at 179.5 fim 
is blended with CH^; the 808 — »7i7 line in Fig. 3f is blended with 
highly excited NH3 (see Fig. 10a); the 524 — > 4i3 (Fig. 3i) and the 
221 ^2i2 lines (Fig. 4q) may have some contribution from CH; 
finally, the spectral feature at 104 yum associated in Fig. 4a with 
the 615 — > 606 line is probably contaminated by OH^. Other con- 
taminated lines marked with dashed-green arrows in Fig. 2 are 
not shown. The OH lines in Fig. 6 have relatively weak contam- 
ination by NH3 in the red A-component of the 84 fim doublet 
(panel e, see Fig. 10c), and by H2O in the blue A-component of 
the 65 fim doublet (panel f). 

The detections summarized in Fig. 2 indicate extreme ex- 
citation in the nuclei of both galaxies. In Arp 220, lines with 
lower level energy up to isiower ~ 600 K in both H2O {7qj — > 616 
and 7i7— »6o6. Fig. 3j, the latter contaminated by excited NH3) 
and OH (ni/29/2 ^ 7/2, Fig. 6i) are detected. NGC 4418 
shows even higher-lying lines, with detected H2O lines as high 
as 818— >7o7 and 7i6— >625 (Fig. 3f,g; /iiowei- Z 800 K). These are 
the highest-lying H2O lines detected in an extragalactic source 
to date. 

The relative strengths of the H2O and OH Unes in NGC 4418 
and Arp 220 show a striking dependence on /iiower- The lowest- 
lying (flower ^ 200 K) H2O lines, especially the 2i2— > loi tran- 
sition at 179.5 yum (Fig. 4q), are stronger in Arp 220 than in 
NGC 4418, but higher-lying fines ai-e stronger in NGC 4418 (see 
Fig. 13). Tfie effect is even more clearly seen in OH (Fig. 14), 
wfiere tfie ground-state lines (upper row in Fig. 6) are mucfi 
stronger in Arp 220, tfie intermediate 113/25/2 — » 3/2 transitions 
at 84 jum (Fig. 6e) fiave similar strengtfis, and fiigfier excited OH 
lines (Fig. 6f,g,fi,i) are stronger in NGC 4418. Tfiis is partially 
a consequence of tfie Cextended component in Arp 220: since Weq 
are calculated relative to tfie observed (nuclear-i-Cextended) con- 
tinuum, tfie fiigfi-lying lines arising from tfie nuclei are reduced 
due to tfie contribution by Cextended to tfie continuum, wfiicfi is 
more prominent in Arp 220. In otfier words, tfie fiigfi-lying lines 
are more continuum-diluted in Arp 220 tfian in tfie relatively 
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Fig. 5. Marginally detected H2O lines in NGC 4418 (upper spec- 
tra in each panel); the corresponding lines from Arp 220 are also 
shown (lower spectra). Model predictions for these lines (for Cco,c in 
NGC 4418 and for C„esi in Arp 220) are also shown (red curves). 



"naked" NGC 4418, wfiicfi is consistent witfi tfie continuum 
models sfiown in §3. On tfie otfier fiand, Cextended in Arp 220 pro- 
duces strong absorption in tfie low-lying lines, tfius increasing 
tfieir Weq relative to tfiat in NGC 4418. 

Nevertfieless, dilution in tfie continuum cannot explain tfie 
fiigfier excitation indicated by tfie very fiigfi-lying lines of H2O 
and OH in NGC 4418. Tfiis is best seen in Fig. 13gfi, wfiere tfie 
line fluxes are normalized relative to tfiat of tfie excited 330 — » 221 
H2O line at 66.4 yum. Tfie latter line (Fig. 3g) is well detected 
in botfi sources and, witfi Slower - 160 K, is not expected to 
be significantly contaminated by absorption of extended low- 
excitation H2O in Arp 220. Lines witfi Eiower above 500 K 
are sfiown to produce, relative to 330^221, more absorption in 
NGC 4418 tfian in Arp 220. Also, as mentioned above, tfie H2O 
8i8^7o7 and 715^625 transitions (Fig. 3f,g) are only detected 
in NGC 4418 (tfie 8o8^7i7 line in Fig. 3f is contaminated by 
NH3, and tfie feature is not detected in Arp 220). Two more para- 
H2O transitions detected only in NGC 4418 are tfie 533— >422 
and 440 — > 431 lines (Fig. 3a,r); tfie lower levels of tfiese lines are 
non-backbone and tfius require fiigfi columns to be significantly 
populated. Otfier H2O lines detected only in NGC 4418 are tfie 
625 — >6i6 and 441 — »432 at 94.7 yum, but tfie Arp 220 spectrum 
at tfiis wavelengtfi sfiows a very broad feature, possibly due to 
pointing or instrumental effects (Fig. 3q). 

Tfie spectrum of NGC 4418 clearly sfiows two intrinsically 
weak H2O lines witfi AK+i - +3, tfie 431 — >4o4 line at 61.8 yum 
and tfie 330— »3()3 line at 67.3 yum (Fig. 3e,fi). All otfier detected 
H2O lines fiave fiigfier transition probabilities and, witfi tfie ex- 
ception of tfie very fiigfi-lying lines, are expected to be strongly 
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saturated. While the ortho-H20 67.3 yum line is also detected 
in Arp 220 within a wing-like feature that probably has some 
contribution from Hj^O 330— > 221, the para-H20 61.8 yim line is 
not detected in Arp 220. The high columns in NGC 4418 are 
also confirmed by comparing the absorption in the ortho/para 
pair44i— »432 and 440 —> 431 (Fig. 3q,r), which have similar lower 
level energies {k, 550 K), wavelengths, and transition probabil- 



ities (A » 0.15 s"'). For an ortho-to-para ratio of 3, one would 
expect the same ratio of 3 for in the optically thin regime, 
but this is as low as 1.15 indicating that both lines are still well 
saturated. 

Even higher-excitation lines of H2O, with Slower ~ 900 - 
1100 K, may be present in the spectrum of NGC 4418 (blue- 
dashed arrows in Fig 2). Among these features, shown in Fig 5, 
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the most likely detection is the 726 ^615 hne in Fig 5c, all 
other features being marginal. In particular, the 92.8 yum feature 
in (f) has a non-Gaussian shape. Nevertheless, all features to- 
gether appear to indicate more extreme excitation in the nucleus 
of NGC 4418. This is plausible, as HCN rotational Hnes with 
Slower ^ 1075 K are detected in NGC 4418 (see below). 

As discussed earlier, dilution in the continuum also lowers 
the Wsq of the high-lying OH lines in Arp 220 relative to the 
values in NGC 4418; however, the most excited OH transitions, 
the ni/2 9/2 ^7/2 and Hi/a 11/2^9/2 lines at 56 and 53 yum, 
respectively, are intrinsically stronger in NGC 4418 according to 
our decomposition of the continuum. Another effect of Cextended 
is also worth noting: the only OH transition observed in emission 
above the continuum, the ni/2 3/2 — > 1/2 doublet at 163 yum, 
is excited through absorption of continuum photons in the 35 
and 53.3 fj.m doublets followed by a cascade down to the ground 
OH state (see App. II by Genzel et al., 1985), and is therefore 
expected to arise in the Cextended component (G-A04). The strong 
OH 163 yum emission in Arp 220 reflects the massive envelope 
around its nuclei. In NGC 4418, the doublet is weaker but still in 
emission, also hinting at the existence of a Cextended component 
around its nucleus. 



2.4. Hl^Oand^^OH 

The spectra around relevant H^^^O and '^OH transitions in 
NGC 4418 and Arp 220 are displayed in Figs. 7 and 8. Close 
spectral features due to H30^, NH2, C3, and CH^ are also in- 
dicated. Hj^O is detected in both sources, though in relatively 
low-lying lines. The highest-lying detected line is the 432— >32i 
transition with /siower ~ 270 K in Fig. 7a, as detection of the 
616 5o5 transition in panel e is rather marginal. The most strik- 
ing feature in Fig. 7 is the relative amount of absorption observed 
in NGC 4418 and Arp 220. While the absorption in the main 
isotopologues is stronger in NGC 4418 than in Arp 220 (Figs. 3 
and 4), the opposite behavior is generally found for the absorp- 
tion in the rare isotopologues. This is best seen in the 331 ^220, 
321 ^2i2, 221 lio (showing redshifted emission in Arp 220), 
3i3 — >2o2, and 221 — »2i2 lines. There are, however, two spectral 
features nearly coincident with the 423 — > 3 12 and 4i4 — > 3o3 Hl^^O 
lines at 79.5 and 114.3 yum (panels d and h, respectively), with 
absorption in NGC 4418 apparently stronger than in Arp 220, 
but the latter is subject to an uncertain baseline. 

The '^OH lines in Fig. 8, specificafly the "*OH 85 and 
120 yum doublets, confirm the '^O enhancement in Arp 220. 
The blueshifted A-component of the '**OH 120 yum doublet is 
strongly contaminated by CH^ 7 = 3 — » 2 in both sources, 
but the uncontaminated redshifted A-component, undetected in 
NGC 4418, is strong in Arp 220. It is also worth noting that 
the two ''^OH 85 yum A-components show rather different ab- 
sorption depths in Arp 220. This A-asymmetry is not expected 
to arise from radiative pumping or opacity effects, because the 
A-components of any transition have the same radiative transi- 
tion probabilities, and the corresponding lower levels tend to be 
equally populated. Indeed, the absorption in the corresponding 
A-components of the main isotopologue at 84 fj.m are very sim- 
ilar (Fig. 6e). The "*OH 85 yum A-asymmetry in Arp 220 is fur- 
ther discussed in §3.2.3. The '^OH A-component at 65.54 jum in 
Fig. 8a is partially blended with two lines of NH2 and, together 
with the expected contribution by H20^ at 65.61 yum, form a 
broad feature with uncertain baseline. The red A-component at 
65.69 fim appears to be stronger in Arp 220 than in NGC 4418. 



11 1 , 11 1 11 , 


331-230 










■ , , 1 , , i 1 , , , 1 " 


"1 , 1 , 1 " 





59.2 59,4 59,6 67,6 67,8 68 75,9 76,2 



d 433- 


1 ' 1 
On -JU 


1 1 1 1 


f 








1 , : 1 , r 


1 , , ■ 1 


" , 1 , : 1 , 1 " 


79.4 


79,6 79,8 82,2 82,5 90,8 91 91, 


1 1 1 1 1 


1 1 1 

-lio 


1 1 1 1 1 1 


_ 1 1 1 1 1 1 1 1 










■ , , ^, 1 , 


" 


- , 1 , , i , , 1 - 





109,2 109,6 114 114,3 114,6 139,5 139,8 140,1 





1 1 1 1 1 _ 

"^12 _ 


-k 




1 1 1 1 L 

n J- 












" 1 , , , 


1 , , , 1 " 


" , 1 , , 


: 1 , , , " 



174 174,4 174,8 183,2 183.6 
Kest (Mm) 



Fig. 7. Spectra around the wavelengths of relevant H^^O lines in NGC 
4418 (upper profiles) and Arp 220 (lower profiles). Model predictions 
for NGC 4418 (C,o,e) with HjO/Hf O = 500, and for Arp 220 (C„es,) 
with H20/H^^0 = 70 are shown in red. 



2.5. HCN 

HCN is a key molecule widely studied in both Galactic and ex- 
tragalactic sources. In NGC 4418 and Arp 220, the low-lying 
lines at millimeter wavelengths have been observed by Aalto et 
al. (2007 a,b) and Wiedner et al. (2004). The V2 = 1 ' vibration- 
rotation band at 14 yum was detected in both sources by Lahuis 
et al. (2007), and rotational emission from the upper vibrational 
state in NGC 4418 has been detected by Sakamoto et al. (2010). 
The V2 = l' 1-type HCN lines from J = 4 - 6 at centime- 
ter wavelengths were detected in absorption toward Arp 220 
(Salter et al., 2008). Among the 15 extragalactic sources de- 
tected in the HCN band, the Lahuis et al. (2007) analysis in- 
dicated that NGC 4418 has the second highest HCN column. 
Recently, Rangwala et al. (2011) have reported the detection of 
several emission/absorption HCN lines with Herschel/SPIRE in 
Arp 220. 

We report in Fig. 9 the PACS detection of HCN absorp- 
tion in pure rotational transitions at far-IR wavelengths in both 
NGC 4418 and Arp 220. In NGC 4418, ah transitions from 
7 = 18 ^ 17 (£io„ei- = 650 K) to 7 = 23 ^ 22 (£io„ei- = 1075 
K) are clearly identified. Furthermore, there are hints of absorp- 
tion up to the 7 = 25 — > 24 transition (fiiowei - 1275 K). In 
Arp 220, the 7 = 20 — » 19 line is clearly contaminated by H30^ 
and NH2, and the 7 = 23 — > 22 one by C3 and possibly also by 
NH2. The 7 = 19 —> 18 line is weak, suggesting that this tran- 
sition is tracing the tail of the Spectral Line Energy Distribution 
(SLED). However, a clear feature is found at the wavelength of 
the 7 = 21 — > 20 line, and the 7 = 22 — > 21 line could show 
some marginal absorption as well. Nevertheless, we consider the 
identification of the 7 = 21 — > 20 line as questionable given the 
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weakness of the nearby rotational HCN lines and the possible 
contamination by other species. 

Fluxes derived from Gaussian fits are shown in Fig. 15 as a 
function of fiiower. where values derived from Herschel/Spire ob- 
servations by Rangwala et al. (2011) are included for Arp 220. 
The HCN excitation in NGC 4418 is extreme, with the SLED 
apparently peaking at 7 = 21 ^20 (/siower ~ 900 K). The 
striking characteristic of the far-IR HCN lines is that they are 
detected in absorption against the far-IR continuum, indicating 
that the dust temperature and far-IR continuum opacities behind 
the observed HCN are high; otherwise the lines would be de- 
tected in emission. We argue in §3 that the H2O, OH, and HCN 
lines in NGC 4418 are tracing a high luminosity, compact nu- 
clear source, denoted as the nuclear core (Ccore)- In Arp 220, the 
HCN SLED as seen by Herschel/Spire (Rangwala et al., 2011) 
peaks aty = 16 — > 15 (Slower ~ 500 K) indicating, like H2O, 
more moderate excitation, but a second, higher-excitation com- 
ponent could be present if the 161.3 yum spectral feature is due 
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Fig. 12. Comparison between the line shapes of the HCO^ 7 = 3 — > 2 
(green histograms) and 7 = 4 — > 3 (blue histograms) transitions toward 
(a) the western and (b) eastern nucleus of Arp 220 (from Sakamoto 
et al., 2009), with selected lines detected with Herschel/PACS (c-i). 
The dashed vertical line indicates our adopted systemic velocity (§2.2). 
The velocity scale for the OH 79 (163) fim doublet is relative to 
the red (blue) A-component, and the systemic velocity for the other 
A-component is also indicated. The velocity of the H'^^O 3o3 — » 2i2 line 
in (c) is also indicated. 



to HCN 7 = 21^ 20. As shown in §3.2.4 and 3.4.1, high HCN 
columns are required to explain the data. 

2.6. NH3 

Ammonia is another N-bearing species widely observed in 
Galactic sources mostly through the pure-inversion transitions 
from metastable (/ = K) levels at centimeter wavelengths 
(see Ho & Townes, 1983, for a review). The rotation-inversion 
transitions lie at far-IR wavelengths, and many of them have 
been detected with ISO toward Sgr B2 (Ceccarelli et al., 2002; 
Polehampton et al., 2007) indicating a hot, low dense molec- 
ular layer interpreted in terms of shock conditions (Ceccarelli 
et al., 2002). Detections of far-IR NH3 lines in extragalactic 
sources were presented, together with a preliminary analysis, by 
G-A04 in Arp 220. In Fig. 10 we show the Herschel/R\CS de- 
tection of high-excitation, far-IR NH3 lines in both NGC 4418 
and Arp 220. 

NH3 has two species, ortho-NHs {K = 3n) and para-NH3 
(K + 3n), with an expected ortho-to-para abundance ratio of 



unity (Umemoto et al., 1999). All rotational {J,K) levels, ex- 
cept those of the K - ^ ladder, are split into 2 inversion dou- 
bling sublevels (a for asymmetric and s for symmetric), which 
for a given J have decreasing level energy with increasing K. 
Radiative transitions are only allowed within /T-ladders, thus 
only collisions can populate the A' > 1 ladders and the absorp- 
tion from different /T-ladders is then sensitive to the gas temper- 
ature. Within a given /T-ladder, both collisions and absorption 
of far-IR photons can pump the J > K non-metastable levels. 
Since the a ^ s (s ^ a) J + \ — >7 transitions of different 
/T-ladders have similar wavelengths, the lines are crowded in 
wavelength and overlapping in some cases, with severe blending 
between the K - Q and K - \ ladders. The RACS domain cov- 
ers the range from the {^,K)a — > {2,K)s (^ 170 jum. Fig. lOg) 
and the {'i,K)s {2,K)a (^ 165.7 ;um. Fig. lOf) lines, up to 
(\Q,K) (9,/:) at 50-51 jum. 

Figure 10 also indicates (partial) blending with lines of other 
species, as H2O (panels a and b), CH^ (panel b), and H30^ 
(panel f). Most (6,K)a — > (5,K)s lines are strongly blended 
with the OH 84 fim doublet. From the comparison with the mod- 
eling described in §3, the (6, K)s — > (5, K)a and especially the 
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(7, K)a — > (6, K)s lines are so blended that a pseudo-continuum 
is expected to be formed, with the consequent uncertainty in the 
baseline subtraction. As also found for H2O, OH, and HCN, 
NGC 4418 shows in NH3 clear indications of higher excita- 
tion than Arp 220, and the absorption in non-metastable lev- 
els relative to metastable ones is also stronger in NGC 4418 
(e.g., the (4,K)s (3,K)a group in panel d). The ortho-NHa 
(8,6)fl (7,6)s is probably detected in NGC 4418, as the 
blended para-HaO 8o8^7|7 line (see also Fig. 3f) is probably 
less strong than the close ortho-H20 813— »7o7 line, and also 
because the (8,6)i — > (7,6)a transition is detected. There are 
also some hints of absorption in the (8,5)s — > (7,5)fl hne. 
In summary, NH3 lines are detected up to (8, K) — > (7, K) in 
NGC 441 8 (Slower ~ 600 K), and up to (7,6)^ (6, 6) in Arp 220 
(Slower » 400 K). 

2.7. Kinematics 

In Fig. 12, we compare the HCO^ 7 = 3-2 (green histograms) 
and y = 4 - 3 (blue histograms) line profiles toward the west- 
em (a) and eastern (b) nucleus of Arp 220 (from Sakamoto et 
al., 2009), with the line profiles of selected lines detected with 
Herschel/PACS (c-i). The HCO+ hnes toward both nuclei, ob- 
served with high angular (0.3") and spectral (30 kms"^) res- 
olution, exhibit P-Cygni profiles indicative of outflowing gas 
(Sakamoto et al., 2009). Specifically, the redshifted spectral fea- 
tures observed in emission have a large velocity extent of ~ 500 
kms"' from the systemic velocity. This prominent redshifted 
HCO"^ emission has its counterpart in several lines detected with 
Herschel/PACS: the H2O 3o3^2i2 line at 174.6 iim (Fig. 12c), 
the [O i] 63 /im transition (d), and the OH 79 fim doublet (e). The 
latter also shows deeper absorption in the blue A-component 
than in the red one, suggesting the occurrence of redshifted emis- 
sion in the blue A-component as well. The line emission fea- 
tures are expected to be formed in gas located at the back/lateral 
sides (i.e. not in front of the nuclei where absorption of the con- 
tinuum dominates the profile), and will therefore be redshifted, 
as observed, if the gas is outflowing. The velocity of the emission 
features (vlsr ~ 5650 - 6000 km s"^) is significantly higher than 
the redshifted velocities measured for the rotating disks around 
both nuclei (up to ~ 5700 kms Sakamoto et al., 2008), in- 
dicating that little of this emission is associated with the rota- 
tion motions. Thus we conclude that moderate excitation hnes 
of H2O, OH, and O i trace the outflow at redshifted velocities 
detected in HCO^. The asymmetric shape of the OH 163 yum 
emission doublet, which peaks at redshifted velocities, is consis- 
tent with this scenario, as is the detection of redshifted emission 
features in the Herschel/SPIRE spectrum of Arp 220 in several 
species including H2O (Rangwala et al., 201 1). 

The blueshifted velocity extent of the OH 163 and [O i] 
145 jum hnes, observed in emission, reaches ~ 400 km s ' from 
the systemic velocity (Fig. 12h-i), and could also trace outflow- 
ing gas at the most extreme velocities. The full velocity extent 
in these lines (» 1000 kms"') is similar to that of CO (3-2) 
(Sakamoto et al., 2009). The full velocity extent of the high-lying 
H2O lines is significantly lower, » 550 kms"'. This is specif- 
ically illustrated in Fig. 12f-g for the 322 ^2ii and 6i6^5o5 
lines, which lie in the range 90 - 82 fim where PACS has rel- 
atively high velocity resolution (130 - 150 kms"'). In contrast 
with the OH 163 and the [O i] 145 jum lines observed in emis- 
sion, which probably trace spatiaUy extended gas, the formation 
of the high-lying absorption H2O lines is restricted to regions 
with high far-lR radiation densities that are optically thick at 
far-IR wavelengths, obscuring the emitting gas behind the nu- 



clei. The absorption in these H2O lines peaks at around central 
velocities, indicating that the lines mainly trace gas rotating on 
the surface of the nuclei. However, due to uncertainties in the 
velocity correction due to pointing shifts (see §2), and the sim- 
ilarity in the blueshifted velocity extent of the HCO^ and H2O 
lines, we cannot rule out a significant contribution to the absorp- 
tion by outflowing gas with velocities up to ~ 150 kms~'. Nor 
can we rule out spatially extended gas inflow, similar to what is 
inferred and discussed below for NGC 4418, since the peak ab- 
sorption in the OH 79 fim doublet is redshifted by ^ 90 km s"' . 

In NGC 4418, the adopted redshift (z = 0.00705) is derived 
from the excited H2O and OH lines, coinciding within » 15 km 
s~' with both the centroid of the [C n] 158 fim line and with 
the redshift inferred from the stellar absorption at optical wave- 
lengths (§2.2). No indication of outflowing gas in the nuclear re- 
gion is found. The excited H2O and OH lines in NGC 4418 are 
indeed significantly narrower than in Arp 220 (Figs. 13 and 14), 
further indicating the relatively quiescent state of the NGC 4418 
nucleus. On the other hand, the line shapes of the lowest-lying 
H2O, OH, and [O i] hnes show an intringuing behavior, oppo- 
site to the outflow signatures in Arp 220: the absorption fea- 
tures in the ground OH and [O i] 63 //m lines are systematically 
redshifted by ~ 100 kms and the [O i] 63 /im line shows a 
blueshifted spectral component in emission, peaking at =s -55 
kms ', i.e. an inverse P-Cygni profile (Fig. 11a). The lowest- 
lying H2O Une at 179.5 //m is also redshifted, but contamina- 
tion by CH^ (2 — > 1) at 179.594 jim makes the case more 
uncertain. Additional clues about this component come from 
the line shapes of the OH 84 and 163 yum doublets: the blue 
A-component of the 84 ^m doublet (£iower = 121 K) shows an 
asymmetric shape with some redshifted excess (Fig. 6e and lie; 
the red A-component is probably contaminated by excited NH3, 
see Fig. 10c), suggesting that some redshifted gas at ~ 50 km s"' 
is significantly excited. The kinematics impUed by the OH 163 
yum emission doublet are unclear due to the low signal-to-noise 
ratio (SNR), as one of the components is apparently blueshifted, 
whereas the other one peaks at central velocities (Figs. 6j, llf, 
and 14). In summary, the ground-state lines of OH and [O i], 
together with the OH 84 and 163 doublets, indicate the pres- 
ence of a component different from the nuclear one in both ex- 
citation and kinematics, most probably extended in comparison 
with the nucleus but much less massive than the Cextended com- 
ponent of Arp 220, and the velocity shifts suggest that this com- 
ponent may be inflowing onto the quiescent, nuclear region. We 
further investigate this scenario in §3.6. 

3. Models 

3.1. Overview 

As the data presented in previous sections have shown, the far-IR 
spectra of NGC 4418 and Arp 220 are dominated by molecular 
absorption, with emission in only some lines and profiles. These 
rich line spectra and their associated dust emission cannot be de- 
scribed by a single set of ISM parameters, but different lines have 
different excitation requirements and are thus formed in differ- 
ent regions of the galaxies. Our approach for both galaxies is to 
fit these different regions and conditions, even though they are 
not spatially resolved with Herschel, with the smallest possible 
number of parameterized components. 

For NGC 4418, we find fliat we need (see Fig. la) (i) a hot 
component (yellow curve, hereafter Chot) that accounts for the 
mid-lR continuum; (;/) a nuclear core (blue curve, Ccore) that 
provides absorption in the high-lying lines of H2O, OH, HCN, 
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Fig. 13. Velocity shifts (a-b), line widths (c-d), and equivalent widths 
(e-f) of the HjO lines in NGC 4418 (left) and Arp 220 (right). Panels 
g-h show the line fluxes normalized to that of the Ssq— »22i transition 
at 66.4 fim (see §2.3 for details). Error bars are l-cr uncertainties from 
Gaussian fits to the lines. 




Fig. 14. Velocity shifts, line widths, and equivalent widths of the [O i] 
and OH lines in NGC 4418 (black symbols) and Arp 220 (blue sym- 
bols). Error bars are l-cr uncertainties from Gaussian fits to the lines. 



and NH3; (///) a warm component (light-blue, Cwarm), that pro- 
vides absorption in moderately excited lines of H2O and OH as 
well as a significant fraction of the far-IR continuum emission; 
and (iv) an extended component (solid green, Cextended), which 
accounts for the low-lying redshifted lines of OH and O i. 

For Arp 220, we need (Fig. lb) (/) a single^ hot compo- 
nent (yellow curve, Chot) that accounts for the mid-IR contin- 
uum emission; (//) the western and eastern nuclear components 
(Cwest and Ceast in blue and light blue, respectively), where the 
high-lying molecular lines are formed; (Hi) the extended com- 
ponent (Cextended, green), which provides a significant fraction 
of the far-IR continuum emission and which is likely associated 
with moderate and low-excitation lines of H2O and OH; and (iv) 
an additional absorbing "halo" component (Chaio), with no asso- 
ciated intrinsic continuum but located in front of the nuclei, is 
required to fit the absorption in the ground-state lines of H2O, 
OH, and O i. 

We use single dust temperatures for every component listed 
above, except for the Cextended component in both galaxies where 
the dust temperature profile is calculated from the balance 
of heating (by the inner components) and cooling (Gonzalez- 
Alfonso & Cernicharo, 1999). In our models, each single dust 

^ Even if the mid-lR arises from the two nuclei, the continuum is 
simulated with just a single component with an effective diameter as 
given in Table 1 . 



temperature component is attenuated by a foreground, screen- 
like shell, which is parameterized by its dust opacity at 25 jum, 
T"25.fer- These screen-like shells are responsible for the silicate ab- 
sorption features imprinted on the various modeled components 
in Fig. 1. The dust temperature profile in Cextended is calculated 
by assuming spherical symmetry and a single nuclear heating 
component located at the center of the modeled source. 

Parameters of the continuum models, and the inferred molec- 
ular parameters, are listed in Tables 1 and 2, respectively. A 
sketch of the modeled sources, showing the approximate spa- 
tial scales of the different far-IR components (excluding the Chot 
component of both galaxies), is shown in Fig. 16. Each compo- 
nent, however, can be interpreted in terms of a single source, as 
implicitly assumed, or alternatively applied to each one of an en- 
semble of Nc smaller clouds of radius Rc that do not spatially 
overlap along the line of sight. The scaling between these two 
approaches is discussed in G-A04. 

The models for Arp 220 are discussed in sections 3.2 (high- 
lying lines) and 3.3 (low-lying lines), while those for NGC 4418 
are developed in sections 3.4 (high-lying lines), 3.5 (mid- 
excitation lines), and 3.6 (low-lying lines). The continuum of 
NGC 4418 is further discussed in 3.7. 
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Fig. 15. HCN line fluxes in NGC 4418 (a) and Arp 220 (b). In Arp 220, 
the fluxes of lines with Eiowei up to 570 K have been taken from SPIRE- 
FTS data (Rangwala et al., 201 1). The curves and triangles show model 
predictions for the Ccorc component of NGC 4418 and the C„est compo- 
nent of Arp 220. The dotted curve in (b) shows results obtained when 
only colUsional excitation is included in the model. The solid curve in 
(a) shows the best fit model obtained with T200 ~ 1; while the dashed 
and dotted curves correspond to T200 = 4 with and without gas-dust 
mixing, respectively. 



3.2. The high-lying lines in Arp 220: the nuclear region fCwest 
and Ceastj 

3.2.1. H2O 

Our models for Arp 220 are similar to those described in G- 
A04, and generated on the basis of the interplay between the 
continuum emission and the molecular line absorption. Since the 
molecular excitation in Arp 220 is high and coUisional excitation 
alone cannot account for it, the excitation mechanism is expected 
to be dominated by absorption of dust-emitted photons in the 
nuclear region of the galaxy (G-A04, Gonzalez-Alfonso et al., 
2008, 2010). The molecular excitation is thus a function of the 
dust temperature, Tj, and of A'h2o/'^v', and in combination with 
the observed continuum provides clues about the general proper- 
ties of the far-IR continuum source and its asociated chemistry. 
If the observed H2O excitation cannot be reproduced by assum- 
ing that the dust in the nuclear component has a temperature of 
"only" Td - K (derived from the observed 25 and 1300 /im 
emission), foreground extinction in the far-IR is included to at- 
tenuate also the 25 foa emission from the nuclei in such a way 
that the far-IR is reproduced with > 85 K. This is Scenario 
2, S2, in G-A04, based on the continuum models by Soifer et 
al. (1999). Evidence for such foreground extinction in the mid- 




-100 100 -500 500 

Scale (pc) Scale (pc) 

Fig. 16. Schematic representation of the modeled sources, showing the 
approximate spatial scales of the different far-IR components (see Fig. 1 
and Tables 1 and 2). The Chot components in both galaxies (that account 
for the mid-IR spectra) are most probably associated with the nuclear 
regions and are not included here. The Chaio component in Arp 220 has 
no associated continuum, and is responsible for the absorption of the 
nuclear continuum (Ceast and C„est) in the ground-state lines of H2O, 
OH, and O i. The plot is an oversimplification of the actual models, 
where the different components are modeled separately to account ap- 
proximately for non-spherical symmetry. 



IR comes from the strong silicate absorption at 9.7 and 18 /jm, 
in the millimeter by self- absorption in CO (2^1) (Downes & 
Eckart, 2007), and in the far-IR from the [O i] 63 yum line that is 
observed in absorption, in contrast to most extragalactic sources 
in which it is observed in emission. 

In the present models, we have simulated the dust emission 
by using a mixture of silicate and amorphous carbon grains with 
optical constants from Draine (1985) and Preibisch et al. (1993); 
the mass-absorption coefficient is 550, 150, and 12.3 cm^/g of 
dust at 25, 50, and 200 jum, respectively, and the spectral index 
is yS « 2 in the far-IR. We have also attempted to disentangle 
the emission from the eastern and western nuclei in Arp 220. 
We have modeled the eastern nucleus (Ceast) as a sphere of di- 
ameter d^ast — 

0.31" (Sakamoto et al., 2008, hereafter Sa08), 
radial opacity (i.e. the opacity along a radial path) at 200 
'Tioo = 2.7, and T^^east = 87 K. The emission from Ceast is atten- 
uated by foreground dust with T25,fgr = 1, yielding 2.5 - 0.02 Jy 
at 25 - 1300 fim -consistent with measurements by Soifer et al. 
(1999) and Downes & Eckart (2007). The corresponding SED 
of Ceast is shown in Fig. lb in light-blue, and accounts for a lu- 
minosity of Least = 3 X 10'^ Lq. For the western nucleus (Cwest) 
we also use as a first approach the size as derived by Sa08 in the 
"disk" approximation, i.e. a sphere with diameter c/west = 80 pc 
(0.23" at 72 Mpc), but allow it to vary to match the observed 
continuum, and leave as free parameters the dust temperature 

west 

and foreground extinction T25,fgr- 
Calculations for the lines were carried out in spherical 
symmetry using the code described in Gonzalez-Alfonso & 
Cemicharo (1999). The models indicate that most high-lying ab- 
sorption lines are formed in the outermost shell surrounding the 
far-IR continuum source. This is illustrated in Fig. 17, where 
the absorbing flux of several lines of H2O, OH, Hj^^O, NH3, and 
HCN, is plotted as a function of the depth of the shell where the 
molecules are located, which is parameterized in terms of the 
continuum opacity at 50 fim (T50) measured from the surface of 
the far-IR source. Since the molecular abundances relative to the 
density of dust are uniform in these models, the column densities 
are proportional to T50. Results indicate that the H2O and OH ab- 
sorption is produced in a thin shell with T50 ~ 1, as extinction, 
thermalization by dust emission, as well as line opacity effects 
make results insensitive to the presence of molecules deeper into 
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Fig. 17. Modeled absorbing flux in Arp 220 of several a) H2O, b) OH, 
c) O, d) NH3, and e) HCN lines as a function of the depth of the 
molecular shell measured from the surface of the continuum source, as 
parameterized by the dust opacity at 50 fim (tso)- Numbers in parenthe- 
sis indicate rounded wavelengths in yum. The total column of molecules 
is proportional to T50, and we have adopted A'h20 = 2 x 10'* x T50 
cm-', and abundance ratios of OH/HjO = 0.5, HjO/Hf O = 100, 
HCN/H2O = 0.3, and HCN/NH3 = 6. The dust temperature is 110 
K. 



the far-IR source. The optically thinner H^^O lines are mostly 
formed in the same shell, though the lines at long wavelengths 
(e.g. the 221 — >2i2 at 183 //m) still have significant contribution 
(< 25%) from the 1 < T50 < 2 region. The case for NH3 and 
HCN is somewhat different, as some lines probe deeper regions. 
The high-lying lines of NH3 lie at A < 100 /im and are also 
formed in the T50 ~ 1 outermost shell, but the low-lying lines 



at A > 100 fim have significant contribution from 1 < T50 < 3. 
The high-excitation HCN lines sample regions deeper than the 
other species (T50 < 4) as their excitation is sensitive to the far- 
IR radiation density at A> 170 yum (see §3.2.4). Based on these 
results, the models for the nuclei shown below use a "screen" 
approach, i.e. the molecules are mixed with the dust but located 
within a thin shell surrounding the nuclei extended up to a depth 
of T5() - 2 for H2O, OH, and their '^O isotopologues, and up to 
T"5() = 4 for NH3 and HCN. The values of the column densities 
of all species below, however, are given for the T50 - 1 outer- 
most shell that dominates the absorption of H2O and OH, and 
are lower limits to the true columns through the nuclear sources. 
In order to estimate the abundances xx relative to H nuclei, we 
have normalized these column densities per unit of T50 to a given 
value of A^H.norm- For a mass-absorption coefficient of 150 cm^/g 
of dust at 50 yum and a gas-to-dust mass ratio of 100, T50 = 1 
corresponds to a column density of A^H.nor 



4 X 10^^ cm 



where H refers to hydrogen nuclei in both atomic and molecular 
forms. This value of A^H.norm has been applied to all species ob- 
served in the far-IR, and define the far-IR "photosphere" where 
most molecular absorption is produced. 

Line broadening is caused by microturbulence with Vtur - 
60 kms"^ and a velocity shift of 130 kms ' through the ab- 
sorbing shell. The latter simulates either the presence of out- 
flowing gas or, more generally, gas velocity gradients across 
the nuclear regions, allowing us to nearly match the observed 
linewidths. Rates for collisional excitation are taken from Faure 
et al. (2007); with the adopted «Ht = 5 x 10^ cm"-' and gas tem- 
perature Tg - 150 K (see §3.2.4), collisional excitation has little 
effect on the calculated fluxes. The line models have three free 
parameters: the dust temperature, Tj, the H2O column, Nh^q, 
and the covering factor of the continuum source, fc- 

The bulk of the H2O absorption is produced in a luminous, 
compact region with high H2O columns, which we identify with 
the double nucleus of Arp 220. A number of previous obser- 
vational studies with high angular resolution have revealed that 
the western nucleus is brighter than the eastern one in the near- 
IR (Armus et al., 1995), mid-IR (by a factor of ~ 3 at 25 jum; 
Soifer et al., 1999), and (sub)millimeter wavelengths (Downes 
& Eckart, 2007; Sakamoto et al., 2008). Given that the excita- 
tion of H2O requires a high brightness continuum source, the 
observed high-lying H2O absorption is tentatively attributed to 
Cwest- Nevertheless, the column densities derived below are in- 
dependent of whether Cwest alone is responsible for the absorp- 
tion in the high-lying lines, or Ceast has a significant contribution 
as well, and the properties we derive for Cwest and Ceast can be 
more generally interpreted as shared by both nuclei. 

In the nuclear region where the H2O lines are formed, 
the models shown below indicate that most of these lines are 
strongly saturated; however, there are still some high-lying lines 
that are sensitive to both A^HjO and T^. For these critical fines. 
Figure 18 shows the observed and modeled fluxes relative to 
that in the 330^221 line at 66.4 /im (Fig. 3g). As mentioned 
in §2, the 330 — > 221 is chosen as the normalization line because 
it is detected with a high SNR and, with iiiowei = 160 K, is not 
expected to be sigiuficantly contaminated by absorption of ex- 
tended low-excitation H2O in Arp 220. On the other hand, the 
330— > 221 transition is still low in energy as compared with the 
high-lying lines, so that the line flux ratios are sensitive to both 
NH2O and Td- Furthermore, the line flux ratios do not depend on 
either the size of the continuum source, or on fc. These ratios 
are plotted as a function of for diff"erent "screen-like" H2O 
columns, ranging from Nn.o = 7.5 x lO''' to 6 x lO'** cm , with 
dashed fines indicating the measured values, and dotted lines the 
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estimated upper and lower limits. We also include in the low- 
est panel the upper limit derived for the Sis ^ loi line. For each 
Td, the foreground extinction is determined by imposing that the 
continuum flux density at 25 yum matches the observed 7.5 Jy 
from the western nucleus (Soifer et al., 1999), but this has httle 
effect on the line flux ratios as the far-IR foreground opacities 
are low and, in any case, the considered lines have similar wave- 
lengths. The minimum Td is 90 K, which is the optically thick, 
lower limit value derived for the western nucleus from the 860 
/im submillimeter continuum by Sa08. 
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Fig. 18. Modeled H2O line ratios in Arp 220 as a function of the dust 
temperature in the western nucleus. The dashed lines indicate the ob- 
served values, and the dotted lines the estimated upper and lower lim- 
its. In the lowest panel, the dotted line is the 3cr upper limit for the 
818 - 707/330 - 221 flux ratio. Squares, triangles, circles, and crosses 
(only for rawest = 90 K) show model results for "screen" H2O columns 
of Whjo = 7'.5x 10", 1.5x 10'^ 3x 10'*, and 6x 10'* cm'^, respectively. 

Results of Fig. 1 8 show that the line ratios scale linearly with 
Td. The highest sensitivity to Td is found for the 615 ^524 line, 
since its lower level is non-backbone (Fig. 2). Among the de- 
tected Unes, the line ratios can be almost equally well repro- 
duced with Td - 90 K and the highest A^HjO = 6 x 10'** cm"^, 
or wifli Td = 140 K and /VH2O = 7.5 x 10'^ cm ^. However, 



flie model Td = 140 K and A^HjO = 7.5 X lO'^ cm^^ yields too 
much absorption in the undetected 818 — » 7o7 line, slightly favor- 
ing Td < 130 K. Concerning the continuum models, if we as- 
sume that the observed high-lying H2O absorption is dominated 
by Cwest, and that the photosphere at 25 /urn has the same size 
as the (sub)millimeter source, a reasonable match to the 20 - 35 
fim SED with Td,^xe!it < 1 10 K is not possible; for example, the 
model with Tj ^est = 90 K yields « 4 Jy at 25 //m with no fore- 
ground extinction and the size derived from the submillimeter 
(Sa08), insuflicient to account for the observed 7.5 Jy at 25 jum 
(Soifer et al., 1999). Furthermore, any reasonable fit to the con- 
tinuum requires foreground absorption at 20 - 35 fim, and the 
absolute H2O fluxes also require a continuum source larger than 
the submillimeter one for ra west ^ 120 K. However, the photo- 
sphere at 25 jum (i.e. the nuclear region with T25 > 1) may be 
larger and colder than the submilUmeter source if heated by the 
central core (Sa08), because the optically thinner submillime- 
ter emission samples warm regions that are obscured at 25 /im. 
For Tj^west = 90 K, the required source diameter is d^^st = 150 
pc, yielding L^est = 6.5 x 10" L©, while for Td^west = 130 K, 
<^west = 80 pc and L„est = 8.2 x lO'' Lq. 

In summary, a range of r^,„e.st = 90 - 130 K and a corre- 
sponding "screen" column per unit of T50 of A'hjO - (60 - 7.5) x 
10'' cm"^ are derived, with the colder (< 120 K) sources yield- 
ing diameters above the (sub)milhmeter observed values. The 
corresponding luminosities are in the relatively narrow range 
Lvieat = (6.5 - 8.2) X 10" Lq. The corresponding estimated H2O 
abundances are;ifH20 ~ (15 - 2) x 10"* relative to H nuclei. 

As the reference model shown in Figs. 3-5 for detailed com- 
parison with data, we use the combination Tj ^est = 110 K and 
A^HjO = 2 X 10^^ cm"^, with complete coverage of the source 
(fc = 1). The corresponding continuum is shown with a blue 
curve in Fig. lb, with the diameter of the source (106 pc) in- 
creased by ~ 30% relative to the submilUmeter source, and with 
a foreground opacity T25,fgr = 1.4. The data also suggests the 
presence of a lower excitation component to attain a better fit 
to the 322— >2ii, 221— >lio, and 3o3— >2i2 lines at 90, 108, and 
174 jum, respectively. Tentatively associating this component 
with the eastern nucleus, we derive A'hjO = 3 x lO'^ cm"^ and 
fc = 0.4. Its contribution to the H2O SLED, shown with light- 
blue curves in Figs. 3 and 4, is expected to be significant only 
for low-lying H2O lines. Continuum and hne parameters of the 
models are listed in Tables 1 and 2. 

3.2.2. OH 

The observed high-lying OH doublets in Arp 220, namely the 
n3/2 9/2 7/2 and 11/2 9/2 at 65 and 53 /im, and flie ni/2 
7/2 5/2 and 9/2 7/2 at 71 and 56 jum, are compared in 
Fig. 6 (lower spectra) with the models for the western and east- 
ern nuclei (blue and light-blue curves). CoUisional rates between 
OH and H2 were taken from Offer et al. (1994). As for H2O, flie 
absorption in these lines is dominated by a source with high OH 
column, which we tentatively identify with Cwest- In the model 
of Fig. 6, we have used the same parameters as for H2O, and 
varied the OH column density to match the observed absorption 
in the high-lying Unes. We derive a OH/H2O abundance ratio of 
0.5 - 1 ; higher OH columns overpredict the absorption in the OH 
53 fim doublet (Fig. 6h). The fit to the high-lying OH lines is sat- 
isfactory except for the underprediction of the width of the 03/2 
9/2 — » 7/2 65 /zm doublet, which is probably contaminated by 
relatively weak lines of H2O (gray curve in Fig. 6f), NH2, and 
H20^. Using the same source sizes as for H2O, we find a similar 
OH/H2O -0.5-1 ratio for = 90 - 130 K. An additional 
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extended region is required to fit the low-lying OH lines, which 
is discussed in §3.3. 

3.2.3. Hl^O and i^OH 

Although the detected H^^O lines arise from low-lying levels, 
we find that these absorption features are mainly produced to- 
ward the nuclear region with the highest column densities. Fig. 7 
shows the best fit for the H2^0 Unes using the same reference 
model as for the main isotopologue, i.e. = 110 K. The best 
fit 112*^0 column density is 3 x 10^^ cm"^, implying an Hj^O- 
to-Hj^O ratio as low as ~ 70, and the upper limit for this ratio 
is estimated to be 95. The H^''0-to-H^'*0 ratio slightly de- 
pends on the dust temperature, and a value of » 100 - 150 is 
found for the model with = 90 K, though the overall fit is 
in this case less satisfactory. Thus our data suggest an enhance- 
ment of '^O in Arp 220 relative to the solar value of ~ 5. This 
is confirmed by the models for '^OH shown in Fig. 8, where 
the same '*'OH-to-"^OH ratio of 70 is used with = 110 K. 
The '''OH-to-"^OH ratio would also increase up to 120 - 150 
K for Td = 90 K. The "^OH 65 and 85 jum Unes are nearly 
reproduced, but not the asymmetry in the '^OH 85 jum dou- 
blet (Fig. 8b). Similar intensity asymmetries of x 40% have 
been also observed in most OH A-doublets by Goicoechea et 
al. (2011) toward the Orion bar, where the lines are observed 
in emission. As Goicoechea et al. (2011) indicate, this is proba- 
bly due to asymmetries in the coUisional rates between OH and 
para-H2 (see also Offer et al., 1994). In our absorption lines, the 
asymmetry would arise from the colhsional excitation of '^OH 
in the ground II3/2J = 5/2 — > 3/2 transition, which favors the 
J - 5/2"^ sublevel where the 85.1 fim component arises. Since 
collisions with ortho-H2 quench the asymmetry, the gas is not 
very warm but should be dense enough so that coUisional exci- 
tation with para-H2 competes efficiently with the excitation in- 
duced by the radiation field. The asymmetry is not observed in 
the main isotopologue (Fig. 6e), indicating that either the higher 
column of OH efiiciently mixes the populations of the sublevels, 
or that OH is tracing warmer gas with increased contribution by 
coUisions with ortho-Hi. However, our model with nn, =5x10^ 
cm"^ and Tg = 150 K does not reproduce it, and several tests 
with higher density and lower Tg were also unsuccesful. Further, 
since the ground-state ^'^OH 120 jum doublet is not reproduced 
with the warm component alone, additional absorption by the 
Cextended and Chalo Components described below, with the same 
'^OH-to-'^OH ratio of 70, is included (Fig. 8c). 

3.2.4. HON 

Our models for the western nucleus of Arp 220 can also repro- 
duce the fluxes of the HCN Unes reported by Rangwala et al. 
(2011) and the 7 = 18 ^ 17 and y = 19 ^ 18 lines detected 
with PACS (Figs. 9 and 15). For HCN, both coUisional excita- 
tion and absorption of far-lR photons in the highest-lying lines 
are required to reproduce the observations. At long wavelengths 
(i.e. relatively low- J lines), colhsional excitation in a warm and 
dense region is the primary excitation mechanism. For moderate 
dust opacities, these lines would be detected in emission above 
the continuum. If collisions are able to excite the molecules up 
to the rotational level where the dust emission becomes opti- 
cally thick and the far-IR radiation becomes strong, absorption 
of far-IR photons in high-/ lines continues to excite the molecule 
to higher-lying levels and the corresponding HCN lines are ob- 
served in absorption. Thus the transition from emission to ab- 



sorption Unes is a measure of the dust opacity and temperature 
of this component, and the high-lying HCN Unes detected in ab- 
sorption are tracing high far-IR radiation densities as weU as a 
warm-dense environment with high HCN column densities. 

A grid of models were generated by varying in the range 
10^ - 3 X 10^ cm-3 and Tg between 110 and 250 K. For each 
value of - Tg, both A'hcn and fc were varied to fit the ob- 
served SLED. In all these models, T^ = 1 10 K was adopted. As 
discussed in §3.2.1, the high excitation HCN lines can be formed 
deeper into the dusty, continuum source than Unes from other 
species (Fig. 17), and we extend our HCN shell up to a depth 
of T5() - 4. Rates for HCN collisions with He were taken from 
Dumouchel et al. (2010), and scaled for the reduced mass of the 
HCN-H2 system to obtain the expected colhsional rates with H2. 
Our most plausible results were found for - (4 - 8) x 10^ 
cm-3 and Tg ~ 150 K, for which A^hcn = (6 - 4) x lO" 
cm"^, respectively. A^hcn is twice these values for Tg = 110 
K. Decreasing nn, below 2 x 10^ cm"^ would involve HCN 
columns above 10 cm^^, implying extreme HCN abundances 
above 2.5 x 10"'' that we consider unlikely. The increase of 
above lO'' cm"-' involves lower HCN columns, but the relative 
absorption in the J - 12 —> 11 line becomes somewhat un- 
derpredicted and the value of fc (or the size of the continuum 
source) is above that derived for H2O. Results for Tg = 150 K 
and = 4 x 10^ cm"^, with the resulting A'hcn = 5.6 x lO'^ 
cm~^, are shown in Figs. 9 and 15. The involved far-IR contin- 
uum source has the same size as for H2O and OH. Similar values 
of A^HCN = (3 - 6) X 10'^ cm-2 are obtained for = 130 - 90 K, 
though the sizes of the continuum and HCN sources are in these 
cases higher and lower than for - 110 K, respectively, and 
similar to those found for H2O and OH. 

In our most plausible models, the HCN/H2O ratio is 0.1- 
0.4 (Table 2), and ^hcn = (1 - 2) X 10"^ In Fig. 15, flie cal- 
culated HCN SLED is compared with results for a model that 
neglects the radiative excitation by dust (dotted curve), show- 
ing that the dust has a flattening efl'ect on the predicted SLED. 
Less than 10% of the total column is stored in levels associ- 
ated with the observed lines; the highest populated levels are 
J = 5,6. Nevertheless, owing to the high continuum opacity 
and dust temperature in the nuclear region, as well as to self- 
absorption in the outermost shells of the HCN region, all HCN 
lines above 7 = 6 ^ 5 are predicted in absorption. The emission 
in the / = 6 -> 5 Une (Rangwala et al., 201 1) is not reproduced, 
probably indicating that the line is formed in a colder and more 
extended region associated with relatively weak continuum. On 
the other hand, the model in Fig. 15 accounts for the observed 
absorption up to 7 = 19^ 18 but fails to explain the absorption 
coincident with the 7 = 21 — > 20 transition at 161.3 foa. If this 
feature is due to HCN, a region as warm as that proposed for 
NGC 4418 (§3.4.1) would be involved. 

Our models also include the HCN excitation through the 

pumping of the vibrationally excited V2 = l' state, but we 
find the effect negligible for the ground-vibrational state lines. 
Evidently, this is a consequence of the low T^ used in our mod- 
eUng. Rangwala et al. (201 1) have proposed that the pumping of 
the V2 = 1 ^ state dominates the HCN excitation, requiring a dust 
temperature of T^ > 350 K. If not completely extincted, this hot 
component would show up at mid-lR wavelengths, and there is 
indeed a hot component (Chot) in our continuum modeling with 
T4 = 400 K (yeUow curve in Fig. lb) attenuated by Ay » 60 
mag, which yields an unattenuated flux at 14 fim of « 2 Jy -twice 
as observed. However, Chot has an effective diameter of only ^ 5 
pc, while the observed HCN absorption must be produced over 
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spatial scales of ~ 100 pc because, otherwise, the far-IR HCN 
absorption would be negUgible. Thus, if the HCN molecules are 
exposed to the radiation of Chot, this radiation will be diluted and 
it remains unclear if it can significantly affect the populations. 
As Rangwala et al. (2011) suggest, the alternative is that there 
is a hot component much more prominent, but very extincted, in 
such a way that it does not show up in the observed SED. The oc- 
currence of this component is plausible, but the HCN molecules 
responsible for the absorption are in front of, and thus apparently 
exposed to the non-diluted, unextincted radiation field character- 
ized by a "moderate" Td in the range 90 - 130 K for Cwest, as the 
continuum models indicate. The associated dust is very likely 
optically thick at 14 ixm and thus would obscure the emission 
from the above buried hot component. Our models for the nu- 
clear region are, however, structureless, and more refined models 
are required to clarify the role of the V2 = 1^ pumping for these 
lines. Finally, we note that our derived HCN column is a factor 
of 20 higher than that derived by Lahuis et al. (2007) from the 14 
yum absorption band; as the authors state, however, their analysis 
is based on a covering factor of unity of the mid-IR continuum 
and their columns are then lower limits. 

3.2.5. NH3 

Our model for NH3 is shown with red lines in Fig. 10 and uses 
the same parameters as for HCN: = 4 x 10^ cm"-', Tg - 150 
K, Trf = 1 10 K, and the same source size. The model treats line 
overlaps, and uses the colhsional rates from Danby et al. (1988). 
The derived column density per unit of T50 is A'nhj « 1.5 x lO" 
cmr^, i.e. HCN/NH3 ~ 4 and;^-NH3 ~ 4 x 10"^. The best fit for 
all lines is found for Tg — 150 K, as for Tg - 250 K the lower- 
lying (3, K) — > (2, K) lines would be underestimated. There is 
still some missing flux in these lines at Tg = 150 K, indicat- 
ing possible low-level absorption in a more extended region. 
Within a given /iT-ladder, excitation through absorption of far- 
IR photons dominates over coUisional excitation. The data are 
consistent with an ortho-to-para ratio of 1, indicating a forma- 
tion temperature of at least ~ 30 K (Takano et al., 2000). Similar 
columns are obtained for = 90 and 130 K, with the appropi- 
ate source sizes as derived for H2O, OH, and HCN. Half of the 
derived NH3 column is stored in the metastable (J = K) levels, 
giving T50 = 1 metastable colunms in agreement with the value 
inferred by Ott et al. (201 1) from the inversion lines for T^x = 50 
K(A^NH3 =0.84x 10'^ cm-2). 

3.3. The low-lying lines in Arp 220: the extended and halo 
components (Cextended and ChaioJ 

3.3.1. H2O 

In Arp 220, the absorption in the low-lying 32i ^2i2 (Fig. 3k), 
221^1 10 (Fig- 4b), 3o3->2i2 (Fig. 4p), and especially the 
ground-state 2i2— > loi (Fig. 4q) lines, require a more extended, 
low-excitation region (Cextended, G-A04). This is consistent with 
the continuum models, as the nuclei cannot account for the bulk 
of the far-IR emission between 50 and 300 /im, and an extended 
component heated in part by radiation from the nuclei, but with 
a possible contribution from extended star formation, is invoked 
(green curve in Fig. lb). In our models, Td in Cextended is calcu- 
lated from the balance between heating (by the far-IR emission 
from the nuclei) and radiative coohng, attaining values between 
40 and 90 K. In the model, Cextended extends up to an uncertain di- 
ameter of ^extended - 650 pc, has a column of A^H ~ SxlO^-'cm"^, 
an average density of tin ~ 10^ cm"^, and a mass of M ~ 3x10^ 



Mo (Table 1). Most of its luminosity (=s 7 X lO" Lq) is likely 
reemission, rather than intrinsically generated. 

With NH20 = 3x10^^ cm-2 and XH2O ~ 4x10"^ the H2O line 
absorption due to Cextended is shown with green curves in Figs. 3 
and 4 only for the fines mentioned above. We have used for this 
component a "mixed" approach (i.e. dust and molecules coexist) 
because of the moderate far-IR opacities. We expect reemission 
in the 3o3 — » 2i2 and 322 — > 3 13 lines (G-A04); in the former line 
reemission is indeed observed but redshifted (§2). On the other 
hand, a transition zone between the nuclear region and Cextended 
is traced by the H2O lines in emission at A > 200 detected 
by Rangwala et al. (201 1), as will be discussed in a future paper. 

The depth of the ground-state 2i2^ loi H2O fine at 179.5 
)Um, however, is not quite reproduced even with the absorption by 
Cextended- The fine is probably contaminated by CH^ (2 ^ 1), but 
its contribution is not expected to completely fill the observed 
absorption. This is not surprising, as foreground low-excitation 
molecular absorption of the far-IR emission from the nuclei is 
not included in the calculations of §3.2. Thus we have also in- 
cluded a low-excitation component, Chaio, with no intrinsic con- 
tinuum but located in front of the nuclei, so that it absorbs the 
far-IR emission from the nuclei. This was denoted as the "halo" 
in G-A04, and it is probably a combination of a part of Cextended 
located in the sight line toward the nuclei, together with more 
widespread, foreground, low-excitation gas. The H2O column of 
this component is uncertain due to line saturation, but it is con- 
sistent with A^HjO ~ (5- 10)xl0'^cm-2 (Table 2). The absorption 
in the [O i] 63 /im fine (Fig. 1 la) is also expected to be produced 
in low-excitation gas, which we identify with the Chaio compo- 
nent. From the observed W^q ~ 45 kms"' (Fig. 14) corrected for 
the estimated fraction of the nuclear 63 yum continuum (36%), 
we obtain Nqo ~ 2.5 x 10'^ cm"^. For an oxygen abundance 
ofxo" = 2.5 X lO"'* (Cartledge et al., 2001), Chaio is character- 
ized by A^H ~ 10^^ cm"-^ and T25,fgr ~ 1 . The corresponding H2O 
abundance is;^'H20 ~ 0.5x10"^ relative to H nuclei, similar to the 
value found in Galactic translucent clouds (Plume et al., 2004). 
The inferred A^h ~ 10^^ cm~^ and 25 jum continuum opacity of 
~ 1 are close to the values obtained from the fit to the contin- 
uum for the absorbing layer in front of the nuclei and of the hot 
component (Table 1). While the column density across the Chaio 
component is only accurate within a factor 2, the above compar- 
ison suggests that Chaio is the origin of the foreground extinction 
of the mid-IR nuclear emission. 

Figure 19b shows the predicted- to-observed flux ratios of 
the detected fines for the full reference model including C„est, 
Ceast, and Cextended- Valucs for the contaminated 7i7— »6o6 and 
2i2 — > loi lines are not plotted. Only two lines, the 431 — >422 (147 
/im) and 615^606 (104 yum) show departures larger than 50%, 
but the small predicted-to-observed value of the latter prob- 
ably indicates that the feature is contaminated by OH"^ (§2). 
Among the undetected lines, the absorption in the 533 — » 422 and 
4i3— »4o4 lines at 53.1 and 187.1 yum (Figs. 3a and 4o) are too 
strong in all relevant models. 

3.3.2. OH 

The far-IR spectrum of H2O has only one ground-state line at 
179.5 fim. Since this line is strongly saturated, the species is bet- 
ter suited to describe the region of highest molecular excitation. 
OH has, however, within the PACS wavelength coverage three 
ground-state doublets with different strengths at 119, 79, and 
53.3 yum (Fig. 2), allowing us to better characterize the Cextended 
and Chaio components. This is illustrated in Fig. 6a-c, where the 



17 



Gonzalez- Alfonso et al.: Herschel/PACS spectroscopy of NGC 4418 and Aip 220 



Table 1. Parameters of the continuum models. 



Source 




Diameter 
(pc) 


(K) 


Td 

at 200 ;um 


(cm^^) 


(10** Mo) 


(Lo) 


25,fgr 


4u 

(Lo) 


NGC 44.1 R 


''hot 






fl 0(17 




S V 1 0"' 

O 1\J 




1 ^ 


q n V 10^ 




C'core 


20 


140 - 150 


1-2 


(0.5 - 1) X 10^5 


0.16-0.33 


(6 - 8) X 10'" 


1.5 


3.1 X 10'" 




c 

^.-warm 


31 


110 


0.5 


2.4 X lO^'* 


0.2 


6 X lO'" 


1.2 


2.9 X 10'" 




^'extended 


200 


90-30 


0.03 


(1.0- 1.3)x 10^' 


0.18 


2.4 X 10'° 





2.4 X 10'" 


Arp220 


Chot 


5.2 


400 


0.007 


3.4 X 10^2 


8 X IQ-^ 


1.4 X 10" 


1.1 


2.0 X 10'" 




^west 


150 - 80 


90-130 


6-12 


(3 - 6) X lO^^-e 


55 - 308 


(6.5 - 8.2) X 10" 


1.4 


3.7 X 10" 




Ceast 


108 


87 


2.7 


1.3 X 10" 


13 


3x 10" 


1.0 


2.1 X 10" 




^extended 


650 


90-40 


0.17 


8 X 10^3 


30 


6.9 X 10" 





6.9 X 10" 



" Component: in Arp 220, C„est (Ccast) is the western (east) nucleus. 

Column density of H nuclei, calculated assuming a mass-absorption coefficient of 12.3 cm^/g at 200 //m and a gas-to-dust mass ratio of 100. 

Estimated mass, assuming spherical symmetry. 
^ Unattenuated luminosity of the component. 
^ Foreground opacity at 25 fim. 
' Attenuated luminosity of the component. 

* The estimated mass of Cwcst, based on the flux density at 860 fim reported by Sakamoto et al. (2008) and on an emissivity index of jS = 2 in the 
submillimeter (ks(,o ~ 0.7 cm^ g"' of dust), is higher than the dynamical mass of ~ 10' Mo estimated by several authors (Downes & Solomon, 
1998; Sakamoto et al., 2008). This probably indicates that )3 - 1 is a better approach in the submillimeter (yielding /cg^o ~ 3 cm^ g"' of dust), 
as estimated for some Galactic massive star-forming regions (e.g. Zhang et al., 2007). 



Table 2. Derived H2O column densities and abundances, column density ratios, and covering factors of the continua. 



Source 




"H2O 

(cm-2) 


^"1120 


OH/HzO" 


HCN/HzO" 


HCN/NH^ 


/c(H20)^ 


fciOHf 


/c(HCN,NH3)^ 


NGC 4418 


^core 


(2-6)x 10'» 


(0.5 - 1.5) X 10-' 


0.3 - 0.6 


0.1-0.3 


2-4 


0.85 


1 


0.7 






~ 5 X 10'* 


(0.5 - 2.5) X 10-' 


-0.4 






1 


1 






^extended 


~ 0.2 X 10'* 


~ 2 X 10"^ 


5-10 






0.45 


0.45 




Arp 220 


^west 


(0.8 - 6) X 10'8 


(0.2- 1.5) X 10-5 


0.4-1.0 


0.1-0.4 


3-5 


1 


1 


1 




^east 


3 X 10'^ 


0.8 X 10-6 


~ 1 






0.4 


0.4 






^extended 


3 X 10'« 


4 X 10-8 


~ 3 






1 


1 






^halo 


(5 - 10) X 10" 


5 X 10-8 


~ 10 






0.7 


1 





" Component (see Table 1). In Arp 220, Chaio is an absorbing, widespread component in front of the nuclear region. 

*' Column densities per unit of dust opacity at 50 jum, T50, for Qore and Cwarm in NGC 4418, and for Cwest and Ceast in Arp 220, with an estimated 
column of H nuclei per T50 of 4 x 10 cm . For the ^extended 

component in both sources, mixed values (i.e., dust and molecules are coexistent 

throughout the source; see text for details). 
° Estimated H2O abundance relative to H nuclei. 
'' Column density ratios. 

Covering factor of the continuum (same value for HCN and NH3). 



highly excited OH in the nuclei is expected to produce weak 
absorption in these ground-state lines, thus requiring a more ex- 
tended region where relatively low-excited OH absorbs a signif- 
icant fraction of the far-IR continuum of the galaxy. Two more 
OH transitions also trace Cextended^ the Ylj,/2 7/2 — > 5/2 doublet 
at 84 fj.m (Fig. 6e), and the Hi/a 3/2 ^ 1 /2 doublet at 163 /j^m. 
The former line, with flower - 120 K, is also too strong to be 
reproduced with the model for the nuclei, and should be signif- 
icantly excited in Cextended- The 163 fim doublet, the only one 
observed in emission, is excited in Cextended through absorption 
of continuum photons in the ground-state 53.3 and 35 //m lines, 
followed by cascade down to the ground OH state (Genzel et 
al., 1985), and thus traces large spatial-scale gas illuminated by 
the strong far-IR continuum (G-A04). We indeed expect the 163 
/im doublet to be in absorption toward the nuclei (Fig. 6j), and 
thus Cextended IS required to produce emission in the line. The 
above pumping mechanism also has observational effects on the 
ground-state 79 fiia line: since every 163 fim photon should be 



accompanied, through downward cascade, by a 79 fim line pho- 
ton, the emission in the 163 jim doublet ensures emission in the 
79 yum transition as well. This is indeed observed in Arp 220 
(Figs. 6b and 1 Id), where the 79 jum doublet shows an emission 
feature at the most redshifted velocities of 300 km s"' with re- 
spect to the systemic one; at lower velocities, one would expect 
a combination of emission and absorption with the latter domi- 
nating for suitable geometries and/or sufficiently large amounts 
of OH in front of the far-IR continuum source. With 163 jum ab- 
sorption in front of the nuclei, and emission from more extended 
regions, the 163 fim doublet is expected to be redshifted in out- 
flowing gas, as observed in Arp 220 (Fig. 6j). 

In our specific model for Arp 220, the required OH col- 
umn and abundance in Cextended are A^oh ~ 10^^ cm"-^ and 
XoH ~ 1.3 x 10"', much lower than towards the nuclear re- 
gion and favoring an OH/H2O column ratio of ~ 3. This ratio 
reverses the value found in the nuclear region where OH/H2O 
is ~ 0.5. The Cextended Component then accounts for the 163 
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Fig. 19. Predicted-to-observed H2O flux ratios in NGC 4418 (upper 
panel) and Arp 220 (lower panel) for the models including all com- 
ponents. See §3 for details. 



jum line, though the redshift is not reproduced as no outflow is 
included in this component, and for the absorption at 84 fim. 
Similar to H2O, however, the full absorption in the ground-state 
119, 79, and 53.3 jum doublets is not reproduced; in particular, 
the absorption in the 79 jum line is predicted to be very weak. We 
again find that the missing absorption can be attributed to Chaio. 
requiring a column of A^qh ~ 8x10^^ cm"^. The contribution of 
this component to the absorption in the ground-state 119, 79, and 
53.3 jum doublets is shown with yellow curves in Fig. 6a-c. The 
OH/H2O column density ratio in Chaio is estimated to be ~ 10, 
but the relative OH and H2O columns in the Cextended and Chaio 
components are uncertain due to the difficulty of disentangling 
the two components. 

The ground-state 35 jum doublet of OH (Hi/i — > Hs/i 5/2 
3/2) has been observed with ISO by Skinner et al. (1997); cor- 
recting their continuum flux density at 35 yum (62 Jy) to the lower 
value of 40 Jy indicated by the Spitzer IRS spectrum in Fig. lb, 
the observed absorbing flux in the 35 //m doublet is 1.3x10"'^ 
erg s"^ cm"2. The joint contribution of Cwest, Ceast, Cextended, and 
Chaio predicts a comparable absorbing flux of « I.I x 10"'^ erg 
s"^ cm"^, dominated by Chaio (G-A04). 



3.4. The high-lying lines in NGC 4418: the nuclear core 

(Ccorej 

3.4.1. HON and H2O 

Similar to Arp 220, NGC 4418 shows strong mid-IR silicate 
absorption (Fig. I, see also Spoon et al., 2001) and absorption 
in the [O i] 63 fim line, indicating important mid-IR extinction. 
We have thus followed a similar analysis as for Arp 220, based 
on the combination of the continuum and line absorption. In 



NGC 4418, line broadening is purely modeled as microturbu- 
lence with Vtur = 60 km s"' . 

The extreme excitation of H2O and HCN provides the most 
stringent requirements for the dust temperature (T^), gas tem- 
perature (Tg), density (^Hz), and column densities in the nu- 
clear region. The high-lying lines of these species arise in 
a compact component that we denote as the "nuclear core" 
(Ccore)- Preliminary models for HCN in the Ccore component of 
NGC 4418 are shown in Figs. 9 and 15, with the HCN shell ex- 
tended up to a depth of T50 = 4. We generated several models 
that included/excluded excitation through collisions, absorption 
of far-IR radiation, and pumping through the V2 = 1 ' state, and 
found that HCN is primarily excited through collisions, with the 
far-IR emission significantly modifying the distribution of popu- 
lations in the highest fines. With the SLED apparently peaking at 
7 = 21^ 20, the 

Ccore model requires nu^ ~ 3x10^ cm ^ and Tg 
in the most probable range of 220 - 250 K, significantly higher 
than in Arp 220. Decreasing by a factor of 2 would increase 
Tg to 300 K. The high density and Tg values are also required 
to reduce the otherwise too strong OH absorption in the ground- 
state fines, to account for the NH3 observations described below, 
and to reproduce the far-IR CO emission (discussed in a forth- 
coming paper). In spite of the high collisional excitation, the 
HCN lines are detected in absorption, indicating that the dust 
in Ccore is warm and at least a fraction of Ccore is opticaUy thick 
at far-IR wavelengths. A reasonable fit to the HCN lines is ob- 
tained with Td = 140 - 150 K and T200 S 1 for the continuum. 
With T200 ~ 1 (solid curve and triangles in Fig. 15), we expect 
the rotational lines 7upper ^ 14 to be in emission. The strength of 
the emission fines would decrease if T200 ~ 4 (dashed curve), and 
further weaken with a screen model for HCN (dotted curve). The 
inferred HCN column per unit of T50 is A^hcn = (6 - 7.5) x 10'^ 
cm 2, i.e. XHCN = (1.5 - 2) X lO"''. As for Arp 220, we find fliat 
the pumping of the V2 = I ' state has Uttle eff'ect on the excita- 
tion and strength of the detected high-7 ground-vibrational state 
lines, as expected from our moderate value for T^ (Sakamoto et 
al., 2010). 

The minimum diameter of Ccore. d = 16.5 pc, is obtained by 
assuming/c = I. The corresponding luminosity is « 6x10^" Lq, 
about 60% of the entire luminosity of the galaxy. The Ccore com- 
ponent cannot be much bigger or warmer as its luminosity would 
exceed the bolometric luminosity of the galaxy (assuming spher- 
ical symmetry for the emitting radiation). We also attempted to 
fit die HCN SLED wifli lower T^ in flie range 1 10- 140 K. At 1 30 
K, the source size is increased by 18% relative to the 150 K case, 
the Ccore luminosity is decreased by 20%, and the 7 = 25 — > 24 
line is underestimated; thus our modeling favors T^ > 130 K. 
With Td - 150 K, most energy is emitted at wavelengths shorter 
than 40 fim, but the flux density at 150- 190 /urn is stiU 4.5 - 3 Jy, 
and this is the continuum absorbed by the HCN. In Fig. la, the 
SED of Ccore (blue curve), with an adopted d = 20 pc, is atten- 
uated by foreground dust with r25,fgr = 1.5. The corresponding 
covering factor for HCN is fc ~ 0.7. With T200 = 1- 2, the radial 
column through Ccore is A^H ~ (0.5 - 1) x 10^^ cm"^, sufficient 
to strongly attenuate the X-ray emission from an AGN located 
at the center. 

The Ccore model was applied to H2O, and results are shown 
in Figs. 3 and 4 (blue curves). The highest excited H2O lines 
are the least saturated ones, which makes them most sensi- 
tive to the column density, and we derive A^H20 in the range 
(3 - 6) X 10'^ cm"^. The H2O abundance relative to H is as high 
asxii20 = (0.5 - 1.5) X 10"^, and the corresponding HCN/H2O 
ratio in Ccore is 0.1 - 0.3. The H2O models also confirm the very 
warm (Tj = 140 - 150 K) source of far-IR emission, though the 
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dust mixed in the H2O shell has a somewhat lower Td - 1 30 K to 
avoid strong emission in the 624^615 hne at 167 /im (Fig. 4n). 
H2O excitation through colhsions has two main effects: first, it 
increases the excitation of the low-lying lines among the lowest 
^ 8 levels of ortho and para-HiO (over that of pure radiative 
pumping), providing a better fit to the observations; second, it 
also significantly aff'ects some high-lying lines that, although ra- 
diatively pumped, are populated from the reservoir of H2O in 
the low-lying levels. This combination of increased collisional 
excitation in the low-lying levels, coupled to the radiative pump- 
ing to upper levels, afl'ects mainly the Sig— >7o7 and 615— » 524 
transitions, which are strengthened by ~ 40% relative to models 
where collisional excitation is neglected. 

Our best fit model for H2O uses T200 = 1 - 1.5; a higher 
value (e.g. T200 ~ 4), still compatible with HCN, would produce 
too strong absorption in the H2O lines at long wavelengths (> 
150 ;um, specifically in the 3()3— >2i2 transition in Fig. 4p). On 
the other hand, the covering factor for H2O is slightly higher 
than for HCN. These points could indicate that the very excited 
HCN is more concentrated than H2O (and mostly OH) toward 
the warmest region of Ccore- Figures 3-4 and 19 indicate that the 
model for Ccore is satisfactory for a number of H2O lines, and 
Fig. 5 shows that the model with A^HjO = 5 x 10^** cm"^ predicts 
absorption comparable to that observed in most of the marginally 
detected H2O fines. 

3.4.2. OH 

The high-lying OH doublets in the 113/2 ladder at 65 and 53 /^m, 
and in the ni/2 ladder at 71 and 56 jim, are compared in Fig. 6 
(upper spectra) with the model for Ccore (blue curves). In this 
model we have also adopted the same parameters as for H2O, 
and the fit to OH shown in Fig. 6 is obtained with an OH abun- 
dance decreased by factor of 2.5 relative to that of H2O. The 
OH 65 and 71 //m doublets are somewhat underestimated. The 
OH/H2O ratio in Ccore of NGC 4418, ~ 0.4, is similar to that 
found in the nucleus of Arp 220. Like H2O, the model for OH 
also favors T200 ~ 1 to avoid too strong absorption in the 163 
jum doublet, and the covering factor for OH is higher than for 
HCN and H2O (Table 2). Similar to Arp 220, all other OH fines 
have important contribution from a more extended region, and 
are discussed below in §3.5 and 3.6. 

3.4.3. NH3 

The Ccore model for NH3, shown with red fines in Fig. 10, uses 
the same parameters as derived above for HCN: hyi^ = 3 x 10* 
cm-3, Td = 150 K, and Tg = 200 and 250 K, but results for 
these two Tg values were found to be very similar. The NH3 
data are mostly compatible with T200 ~ 1- The derived column is 
A^NH, « 2.5x10" cm-2, i.e. HCN/NH3 « 3 and^NH, ~ exlO"''. 
In comparison with the model for Arp 220, the higher Tj in NGC 
4418 translates into a higher excitation of non-metastable levels 
along a given A'-ladder, and the higher Tg impUes higher exci- 
tation across ^-ladders. The former effect is responsible for the 
asymmetry between the 124.6 and 124.9 yum features in Fig. lOd, 
which is not fully accounted for by the model; nevertheless, the 
effect is sensitive to the overlap between the K = and K = I 
lines. On the other hand, the model accounts for the absorption 
in the (8, 6)s — » (7, 6)a fine and, with the contribution of H2O, 
for the feature at 63.47 fim, though the marginal absorption in 
the (8, 5)s — > (7, 5)fl transition is not quite reproduced. As for 
Arp 220, the i3,K)a — > i2,K)s feature (panel g) is not fully 



reproduced, suggesting a more extended, lower-excitation com- 
ponent for the low-lying fines. 

3.4.4. and ^^OH 

Relative to the main isotopologues, the li\^0 and ^^OH lines in 
NGC 4418 are weaker than in Arp 220. We have modeled these 
lines in NGC 4418 with the same nuclear Ccore model as for H2O 
and OH, and modeled results are compared with observations in 
Figs. 7 and 8. For Hj^O, the model for the Ccore component (red 
curves in Fig. 7) that reproduces the 331 ^ 220 and 32i — > 2i2 lines 
uses a H20-to-H2^0 ratio of 500. While the model predictions 
are consistent with most observed li\^0 features, the relatively 
strong H2**0 423— >3i2 line at 79.5 //m (Fig. 7d) is underesti- 
mated by the model. On the other hand, the "^OH lines are only 
marginally detected in the spectrum of NGC 4418. Our model 
(red lines in Fig. 8) also uses an OH-to-^^OH ratio of 500, with 
the result that the undetected red 120.2 /im A-component and 
the ^^OU 85 jum doublet are overpredicted. Still, the low SNR 
^^OH lines at 65 and 85 fim are affected by an uncertain baseline; 
owing to additional uncertainties in the OH column and distribu- 
tion, we estimate a lower limit for the OH-to-'^OH ratio of 250. 
In summary, the "'0-to-'*'0 ratio in NGC 4418 is estimated to 
be > 250, with a best fit overall value of ~ 500. 

3.5. The mid-excitation H2O lines in NGC 4418: ttie warm 
component (Cwarmj 

Figures 3 and 4 show that the model for Ccore underesti- 
mates the absorption in some mid- and low-lying lines of H2O 

(330— > 221, 321 — >2|2, 322— >2|i, 221 llO, 4|4— >3o3, 3i3— >2o2, 

and 2i2— > loi), and Fig. 1 indicates that the SED of the atten- 
uated Ccore cannot account for the bulk of the far-IR emission. 
Both points indicate the occurrence of another nuclear compo- 
nent that we denote as the warm component, Cwarm- 

The parameters for Cwarm are more uncertain than those for 
Ccore- Assuming that the continuum emission from Cwarm is 
attenuated by the surrounding extended component (Cextended) 
described below, the dust temperature is ~ 85 - 110 K and 
d ~ 30 - 40 pc (Table 1). The model for the continuum and 
H2O lines of Cwarm IS showu with light-blue curves in Figs. 1, 3, 
and 4. The H2O fines mentioned above can be reproduced with 
NH2O ~ 5 X 10^* cm"-^, much lower than toward Ccore- With an 
estimated T200 ~ 0.5, the model for Cwarm produces emission in 
some lines at /I > 130 /im. The H2O abundance is several x 10^^. 

The model for Cwarm is applied to OH in Fig. 6 (light-blue 
curves), with OH/H2O = 0.4. The main effect of Cwarm is to in- 
crease the absorption in the 119 and 84 //m doublets, and to gen- 
erate reemission in the 163 yum doublet. Nevertheless the joint 
contribution of the Ccore and Cwarm components predicts the 163 
l^m doublet in absorption, in contradiction to Fig. 6j, thus in- 
dicating the presence of a more extended component described 
below. 

3.6. The low-lying lines in NGC 4418: the extended 

component (Cextendedj 

Figure 11 compares the [O i] 63 and 145 jum lines and the 119, 
79, 84 and 163 jum OH lines in NGC 4418 and Arp 220, empha- 
sizing the different kinematics found in both sources. In Arp 220, 
the OH 163 fim doublet, and the emission features in the [O i] 
63 and OH 79 fim transitions, are redshifted relative to the sys- 
temic velocity, and thus reveal outflowing gas though affected 
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by foreground absorption. The broad [O i] 145 //m line supports 
this kinematic scenario. In NGC 4418, however, the line pro- 
files are narrower than, and in some aspects opposite to those 
in Arp 220. The main spectral characteristics of these lines in 
NGC 4418, as discussed in §2.7, are: (0 the [O i] 63 jum absorp- 
tion is rec/shifted by 110 kms"', and an emission feature is 
detected at blueshifted velocities; (ii) the absorption in the OH 
119 and 79 jum lines is also redshifted by « 100 kms~\ with 
some indications of higher redshift in the 119 jim than in the 79 
fim doublet (Fig. 14); (;;;) a redshifted asymmetry is seen in the 
red A-component of the 84 /im doublet; (iv) the 163 yum doublet, 
with relatively low SNR, is detected in emission with one of the 
components apparently fo/weshifted with a velocity comparable 
to that of the [O i] 63 fim emission feature, though the other one 
peaks at central velocities. The lowest-lying H2O Une at 179.5 
yum is also redshifted, but it is contaminated by CH"^ (2 -> 1). 

On the other hand, the high-lying H2O and OH lines in 
NGC 4418 do not show indications of any outflow, but peak at 
the rest-frame velocity as derived from for the optical absorp- 
tion/emission lines. Only the ground-state [O i] and OH lines 
are redshifted, together with some indications of a lower red- 
shift in the excited OH 84 /im line. This indicates that this kine- 
matic component is tracing a relatively low-excitation compo- 
nent around the galaxy nucleus. Near-IR observations (Evans 
et al., 2003) do not show evidence for a nearby "contaminat- 
ing" companion where these lines could be formed (the closest 
galaxy, VV 655, is at ~ 25 kpc), but a linear, compact struc- 
ture with fainter structures extending radially out to ~ 1" (140 
pc) in radius is apparent. This morphology may suggest that the 
redshifted spectral features are formed in (the inner region of) 
these extended streamers, thereby suggesting the possibility of 
an envelope inflowing onto the nucleus of NGC 4418. 

We first constrain the column of the redshifted absorbing 
component in front of the nucleus by calculating the colunm of 
O I from the Weq of the 63 fim absorption feature. Assuming 
that the line is optically thin at all velocities, that the shell cov- 
ers the whole far-IR continuum source at 63 jum (43 Jy), and 
ignoring reemission in the Une, we derive A^oi = 5.3 x 10^^ 
cm"^. This is a lower Umit, as models described below indicate 
that there is significant reemission from the envelope (i.e. from 
Cextended) in both the 63 jum continuum and in the line; we esti- 
mate from the models a more accurate A'oi = (1.4 - 1.8) x lO'' 
cm"^. Assuming a gas-phase O i abundance typical of translu- 
cent clouds (xoi = 2.5 X lO"'' relative to H, Cartledge et al., 2001 ; 
Lis et al., 2001), A'h ~ 6.5 x 10^^ cm"^ or Av ~ 30 mag. Since 
oxygen depletion may be enhanced for those columns (Bergin 
et al., 2000), the allowed range is Ay ~ 30 - 60 mag. This ex- 
tinction is at least half of that required to account for the siUcate 
absorption at 10 /xm. 

One important parameter for the inflowing scenario is the 
amount of gas that is inflowing on the lateral and back sides of 
the nucleus. In order to constrain the geometry, calculations in 
spherical symmetry are performed and compared with data. The 
model for this extended component (Cextended in Tables 1 and 
2) simulates a spherical envelope with inner and outer radii of 
dinner ~ 20 and /?outer ~ 100 pc (0.15" and 0.7") inflowing onto 
the warm nucleus. The value of Rima is chosen to match the es- 
timated outer radius of Cwarm- The line profiles are best matched 
with a density profile that varies as ~ r^^, with a highest value 
of Wh = 2.5 X 10^ cm"^ in the neighborhood of the nucleus. At 
^outer, «H ~ 10^ cm"^, and thus the value for the outer radius is 
not critical. While the [O i]63 fim absorption can be reproduced 
by either a turbulent shell inflowing with uniform v - -110 
km s or by assuming a radial velocity gradient across the en- 



velope, the OH shifts in Fig. 14 suggest a progressive gas de- 
celeration with decreasing distance to the galactic nucleus, with 
the thickest 119 fan line, tracing the outermost regions, showing 
the highest Vshitt, and optically thinner (79 yum) or excited (84 
jum) lines peaking closer to the rest velocity. We thus have em- 
pirically adopted a velocity field with -30 kms ' at Thinner and 
dv/dr = -2.1 km s"^ pc~' (i.e. -200 kms~' at/?outer). roughly 
simulating gas in nearly free-fall at the outermost radius, and 
subsequent deceleration due to radiation pressure (see §4.2.3) 
and/or interaction with orbiting gas, with gas approaching the 
systemic velocity in the neighborhood of the nucleus. 

Our model for OH, shown in Fig. 11, includes both the nu- 
clear component (the combined emission of the Ccore and Cwarm 
components is shown with blue lines) and the inflowing com- 
ponent (Cextended. green Unes). An interesting feature of the line 
profiles is that, at central velocities where the excited lines show 
deep absorption (and the OH column densities are high), the ab- 
sorption in the ground 119 and 79 yum lines is relatively weak. 
The model for the nucleus (Ccore +C„arm) predicts too much ab- 
sorption in these Unes even adopting the high density of 3 x 10^ 
cm"^ and Tg = 200 K for Ccore • Nevertheless, a reasonable match 
to the profiles involves partial cancellation between the absorp- 
tion toward the continuum source and emission (through direct 
scattering at 1 19 jum, and pumping through 53 and 35 fim for the 
79 fim transition) from Cextended around the nucleus. This process 
also accounts for the emission in the 163 fim doublet, as every 
163 fim photon should be accompanied, through downward cas- 
cade, by a 79 fim line photon. In the model, the 163 fim com- 
ponents are blueshifted by -45 kms"', but the observed shift 
is unclear as discussed in §2.7. On the other hand, the adopted 
velocity gradient predicts absorption at ~ -50 km s"' in the 84 
fim excited line, nearly matching the observed asymmetric 84 
fim line shape. The model for OH is then consistent with an ex- 
tended envelope and favors the scenario of an inflow extended on 
the lateral sides, though the lack of blueshifted emission in one 
of the A-components of the 163 yum doublet leaves the case in- 
conclusive. The OH column across Cextended is A^OH ~ 1.5 X 10"" 
cm"^, andxoH ~ (1-5 - 3) X lO"'. From the absorption in the 
H2O 2i2— > loi Une in Fig. 4q, which is contaminated by CH"^, 
we estimate an uncertain A^H20 ~ 2 x 10'^ cm"^. 

The model for OH has been applied to O i. The absorption 
component of the 63 yum Une is reproduced with an O VOH abun- 
dance ratio of 10^, close to the value expected for translucent 
clouds. We first ignore the possible contribution by the nucleus 
to the emission feature of the [O i] 63 fim line. The resulting line 
shape is compared with the observed profile in Fig. 11a (green 
curve), and shows that (/) the velocity field and density profile 
adopted for OH roughly match the absorption in the [O i] 63 
fim line, though the modeled absorption is somewhat overesti- 
mated at low velocities; (//) in our spherically symmetric model, 
the emission feature is formed in the inflowing envelope through 
scattering of 63 fim photons, overestimating by 70% the flux of 
the observed spectral feature (325 Jy kms~^ between -200 and 
km s '). Our model for O i thus supports an inflow with wide 
opening angle, but with the amount of gas inflowing on the lat- 
eral sides being lower than in front of it. 

The above solution for the blueshifted emission in the [O i] 
63 fim line, i.e. that the feature is formed in the inflowing enve- 
lope, is, however, not unique, because of the possible contribu- 
tion by the nucleus to this feature. Indeed, the [O i] 145 yum line, 
though weak and too narrow as compared with the instrumental 
spectral resolution, is detected in the PACS spectrum (Fig. 1 lb), 
and the line should be formed in the nuclear region. Here we at- 
tempt to address the conditions under which the emission in the 
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Fig. 20. Modeled [O i] 63 jum emission above the continuum, integrated 
over blucshiftcd velocities (i.e. between -200 and km s"' with respect 
to the systemic velocity), for the nucleus of NGC 4418 as a function of 
the H2 density. The dashed line indicates the observed value. 



[O i] 63 j-im line may be also formed in the nuclear region. We 
have then calculated the expected nuclear emission in the [O i] 
63 /im line as a function of the assumed nuclear gas density. In 
all models, the nuclear O i column is fixed to match the observed 
[O i] 145 /im line. We have used a "mixed" approach as the con- 
tinuum opacity at 145 /im is moderate and the inner regions of 
the nucleus can contribute to the [O i] 145 jum line. We have 
adopted the dust model for the C„arm component with = 110 
K, and the gas at Tg - 150 K. Collisional rates between O i and 
H2 are taken from Jaquet et al. (1992). Results for the blueshifted 
(i.e. integrated between -200 and kms"') emission are com- 
pared with the observed value in Fig. 20. At the relatively low 
density of 10^ cm"^, the excitation of the [O i] 63 yum line is be- 
low Td, and the line is predicted in weak absorption. However, 
the absorption changes to emission for higher densities, and the 
blueshifted integrated flux attains a value comparable to the ob- 
served one at ~ 3 X 10* cm"^. Thus a significant contribution 
by the compact nucleus to the [O i] 63 //m emission feature is 
not ruled out from the point of view of the radiative transfer. For 
MHz =3x10* cm"^, the "mixed" O i column is A^o" ~ 1-5 x 10^" 
cm"^, aiidx(y> ~ 3 x 10"^. 

The case of the [O i] 63 //m emission appears then somewhat 
ambiguous on the basis of these radiative transfer models. Still, 
the OH models and the requirement of an extended region to ac- 
count for the emission in the OH 163 /zm doublet, favor Cextended 
as the main contributor to the [O i] 63 jum emission feature, but 
the relative contribution by the nucleus is hard to constrain. We 
further discuss this issue in §4.4. 

3.7. The continuum of NGC 4418 

Our continuum fit to NGC 4418, shown in Fig la and with pa- 
rameters listed in Table 1, is composed of: (i) a hot component 
(Chot> yellow curve) with Tj = 350 K and an effective diameter 
of d X 5 pc, attenuated by Ay = 70 mag. It accounts for the mid- 
IR emission including the silicate 10 //m absorption. (//) The nu- 
clear core (Ccore. blue curve) with Tj = 150 K and d -20 pc, at- 
tenuated by Ay = 8 1 mag. It is required to explain the high-lying 
molecular absorption observed in H2O, HCN, OH, and NH3. {Hi) 
The warm component (C„ai,n, light-blue curve) with - 110 
K and d - ?>\ pc, attenuated by Ay — 65 mag. The parameters 
of this component are somewhat uncertain and could include re- 
gions with lower ~ 85 K. Cwam mainly accounts for some 
mid-excitation lines of H2O and OH. (iv) The extended com- 
ponent (Cextended, soHd grccu curve) with non-uniform Td and 
d 200 pc, which mainly accounts for the redshifted OH ab- 



sorption and for the emission in the OH 163 fim doublet. The 
emission from the above components stiU underestimates the 
observed continuum at A > 100 yum, and therefore we have 
added a more extended, cold component with d = 650 pc and 
Jrf = 50 - 25 K (dashed-green curve in Fig. la). This cold com- 
ponent, however, has no obvious molecular counterpart at far-IR 
wavelengths and is included in order to obtain a reasonable fit 
to the SED at long wavelengths, though the total emission (red 
curve in Fig. la) still underestimates the (sub)millimeter emis- 
sion. Possible causes are the wavelength dependence of the dust 
absorption coefficient, and the uncertain dust opacity and size of 
the nuclear core at (sub)millimeter wavelengths. 

The unattenuated luminosities of the hot, core, and warm 
components are 8.5 x lO'*' (i.e. 85% of the bolometric lu- 
minosity), ~ 7 X 10'^ Lo, and 6 x 10'° Lq, respectively. The sum 
of the three exceeds the bolometric luminosity, implying that, 
unless there are important geometrical effects that make us un- 
derestimate the true galaxy luminosity, most emission of the core 
and warm components is not intrinsically generated, but largely 
reemission due to heating from inside. Nevertheless, clear depar- 
tures from spherical symmetry are inferred from our modeling, 
as e.g. the extinction alfecting Chot is much lower than the extinc- 
tion through Ccore> so the latter component can not be interpreted 
as a homogeneous shell surrounding Chot- 



4. Discussion 

4.1. Model uncertainties 

Our schematic models adopt uniform densities and uniform dust 
temperatures for the nuclear components, which is clearly a 
rough approach as it ignores their inner structures and dumpi- 
ness. The use of spherical symmetry is also certainly inaccurate, 
especially in NGC 4418 where a linear structure for the nuclear 
region is observed in the near-IR (Evans et al., 2003). Thus our 
models should be considered equivalent spherical/uniform ap- 
proaches to complex structures. 

Toward the nuclear region of Arp 220 and NGC 4418, our 
column density estimates are dust-limited, in the sense that the 
line absorption traces the outermost regions where Td < 1. In 
this far-IR photosphere, we estimate that our derived columns 
per unit of T50 are probably accurate to within 50% for a given 
Ta. Conversion to abundances is more uncertain. We have nor- 
malized the inferred column densities to a total H column of 
A^H.norm = 4 X lO^"* cm"^. Howevcr, the above H column de- 
pends on the gas-to-dust ratio as well as on the assumed prop- 
erties of dust grains. Furthermore, the HCN absorption is ob- 
served at wavelengths longer than most H2O and OH lines, 
and thus traces higher H columns than the latter species. The 
molecular abundances are thus reliable only to within a factor of 
2-3. Nevertheless, abundance ratios between different molec- 
ular species are based on the same approach, and are therefore 
more accurate. 



4.2. Tlie nuclear regions of NGC 4418 and Arp 220 
4.2.1 . A hot core chemistry 

The far-IR Herschel/PACS observations of NGC 4418 and 
Arp 220 reveal extreme excitation in the key molecular species 
H2O, OH, HCN, and NH3 . The lines are observed mostly in ab- 
sorption against the continuum, and our analysis indicates that 
very high column densities and abundances of these molecules 
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surround (and are associated with) compact, far-IR thick, lumi- 
nous continuum components. 

One of the most intringuing results of our analysis is the high 
abundance derived for HCN, (1 - 2) x 10"*' in both Arp 220 and 
NGC 4418. At milUmeter wavelengths, HCN has been found 
to be enhanced over HCO"^ in ULIRGs (Gracia-Carpio et al., 
2006) and in Seyferts over starbursts (Imanishi et al., 2007), 
and this has been proposed to be a diagnostic of AGN activ- 
ity. Costagliola et al. (201 1), however, argue that for high densi- 
ties and columns HCO^ is expected to be stronger than HCN in 
XDRs (Meijerink et al., 2007). Models by Meijerink & Spaans 
(2005) predict HCN abundances and cumulative column densi- 
ties in XDRs that fall by far too short to explain the values we 
derive in both NGC 4418 and Arp 220. Likewise, the peak HCN 
abundances predicted by Lepp & Dalgarno (1996) for XDRs, 
a few X 10"^, are about two orders of magnitude lower than the 
values we derive. 

The high HCN abundances in both sources are very likely 
tracing a hot-core like chemistry, as has been proposed by pre- 
vious authors, i.e. by Lahuis et al. (2007) from the HCN 14 /xm 
absorption in these and other (U)LIRGs, Martin et al. (2011) 
for Arp 220 based on multispecies observations at 1.3mm, and 
Aalto et al. (2007a) and Costagliola & Aalto (2010) to explain 
the luminous HC3N emission in NGC 4418. HCN is well known 
to be enhanced in Galactic hot cores, as in the Orion hot core 
where the abundance is ~ 10"^ (Schilke et al., 1992), about 
two orders of magnitude higher than in cold molecular clouds. 
Time-dependent chemical models for hot cores by Rodgers & 
Charnley (2001) predict a linear increase of the HCN abundance 
with time for Tg - 300 K, attaining values of 10"^ relative to 
H2 after ~ lO*" yr of chemical evolution. The authors describe 
this model result as "pathological" as such high abundances are 
not found in Galactic hot cores, and suggest that these regions do 
not survive for more than a million years. However, extragalactic 
nuclear regions may survive long enough to develop such high 
HCN abundances. 

The proposed evolved hot core chemistry finds support in 
the HCN/NH3 ratio of 2 - 5 estimated in both Arp 220 and 
NGC 4418. After > 10^ yr of chemical evolution, the models 
by Rodgers & Charnley (2001) predict that a similar value of the 
ratio can be reached provided that the gas temperature is high 
enough, though HCN is predicted to continue increasing with 
time. The ratio found by Rodgers & Charnley (2001) at such 
late stages of hot core evolution is independent of whether NH3 
was injected into the gas phase from mantles or formed directly 
in the gas phase. 

The hot-core scenario is also quite consistent with the warm 
dust temperatures derived in the nuclear regions of NGC 4418 
and Arp 220, and with the derived high H2O abundances. In the 
Orion hot core, Melnick et al. (2010) derive H2O abundances 
of ~ 10"^ from Herschel/HIFI observations. Theoretically, an 
"undepleted chemistry" where grain mantles are evaporated is 
expected to increase the gas-phase H2O abundance up to values 
> 10"^ (e.g. Charnley, 1997), i.e. orders of magnitude relative to 
the values inferred in cold molecular clouds and PDRs (Melnick 
& Bergin, 2005), where H2O ice is expected to be efficiently 
formed upon hydrogenation of atomic oxygen depleted on man- 
tles in cold environments (Bergin et al., 2000; Hollenbach et al., 
2009). If the gas attains sufficiently high temperatures (Tg ~ 300 
K), H2O will be further enhanced in the gas phase via neutral- 
neutral reactions (e.g. Ehtzur & de Jong, 1978). 

The current models underpredict the low-/ lines of HCN ob- 
served by Aalto et al. (2007a) in NGC 4418, who derived densi- 
ties of ~ 10^ cm"^; we derive much higher densities, ~ 3 x 10^ 



cm~^, for the high-7 lines. Likewise, our models for the high- 
J fines in Arp 220 underpredict the 7 = 6 ^ 5 measurement 
of Rangwala et al. (201 1). It is plausible that these sets of lines 
probe different components, with the higher-7 lines tracing very 
dense gas close to the regions of high far-IR emission with tem- 
peratures of Ta > 100 K, while the low-/ lines trace a surround- 
ing or interclump component with more moderate densities. No 
lines of HNC are found in the PACS spectra of NGC 4418 or 
Arp 220, while at millimeter wavelengths HNC is overluminous 
relative to HCN (Aalto et al., 2007b). 

The high OH abundances in both NGC 4418 and Arp 220 
do not easily fit into the hot core scenario. Observationally, OH 
observations of Orion BN/KL by Melnick et al. (1990) indicate 
that the emission/absorption is mostly produced by the shocked 
gas associated with the outflow, with little contribution, if any, 
by the hot core. Recently, Wampfler et al. (201 1) derived a very 
low OH/H2O ~ 10"^ abundance ratio in the warm (T > 100 
K), quiescent component of the Galactic high-mass star-forming 
region W3 IRS 5. Models by Charnley (1997) indeed predict a 
very low OH/H2O ratio in hot cores as most oxygen at high T is 
expected to be chaneUed into H2O. 

4.2.2. X- and cosmic rays 

Another component, with chemistry dominated by photopro- 
cesses and/or cosmic rays, is thus required to account for 
the OH enhancement in the nuclear region of both galaxies. 
Multitransition observations of OH in the Orion bar and chem- 
ical modeUng by Goicoechea et al. (2011) reveal a dense/warm 
FDR component with an OH abundance of a few X 10"^. 
However, the OH enhancement is expected to be restricted to 
the external Ay < 1 layers where OH/H2O > 1 due to H2O 
photodissociation by the FUV field (Goicoechea et al., 201 1). In 
NGC 4418 and Arp 220, we require widespread OH to account 
for the observations and, if PDRs alone are invoked, the surfaces 
of the nuclear regions in both sources (with typical sizes of 40 
and 100 pc) would have to be covered with such high-density, 
low-Ay PDRs with a covering factor of > 0.5 and, for typical in- 
dividual PDR columns of A^oh ~ afew x 10^^ cm~^ (Goicoechea 
et al., 201 1), > 100 of such regions are needed along the line of 
sight. The impact of FUV photons on large amounts of gas can 
be enhanced under special geometries, as in the walls of cavities 
excarvated by outflows (Bruderer et al., 2010); OH is boosted in 
outflows where H2O is photodissociated by strong FUV fields 
(Melnick et al., 1990; Spinoglio et al., 2000). However, again a 
large number of overlapping outflows would be required across 
colunms of several X 10^^ cm"-^. Furthermore, the outflow signa- 
tures in Arp 220 are only observed in fines of moderate excita- 
tion, and NGC 4418 does not show any hint of outflowing gas, 
indicating that high-velocity (as compared with Galactic stan- 
dards) outflows are not responsible for the observed chemistry. 
Nevertheless, the presence of Orion-fike outflows with velocities 
of ~ 30 kms"' is not ruled out. 

A promising scenario to explain the high OH columns and 
abundances is through the efl^ect of cosmic/X-rays, which pen- 
etrate in high- Ay regions much more deeply than UV photons. 
The OH/H2O ratio as predicted by XDR and CR models of high 
density molecular regions is > 1, and;t'0H easily exceeds lO"*" 
(Maloney et al., 1996; Meijerink et al., 201 1). Large amounts of 
OH can be formed through dissociative recombination of HsO"^, 
which in turn is formed through the sequence O"^ — » OH"^ — > 
H20'^ HsO"^ (Maloney et al., 1996); as we wiU report in a 
forthcoming paper, these molecular ions are indeed detected in 
absorption from excited levels in both sources. However, both 
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NGC 4418 and Arp 220 have strong X-ray deficits relative to 
the FIR emission (Maiolino et al., 2003; Iwasawa et al., 2005). 
From the X-ray flux reported by MaioUno et al. (2003) from 
NGC 4418, Fx(03 - 8keV) = 1.5 ergs-'cm-^, and assum- 
ing a point-like X-ray source located at the center of our Ccore 
component, the X-ray flux at the surface of Ccore ('' ~ 10 pc) 
is ~ 0.1 ergs"' cm"^; a similar estimate for Arp 220 with the 
fluxes reported by Iwasawa et al. (2001) yields an X-ray (0.5 - 10 
keV) flux at the surface of the western nucleus (r ~ 50 pc) 
of ~ 0.5 ergs ' cm"^. Since the lines observed with PACS are 
formed around far-IR thick components with Nn > 5 X 10^"* 
cm"^, any putative strong and buried X-ray source would be 
severely extincted, and the above distance-corrected estimates 
are probably lower limits due to foreground absorption but still 
relevant within a factor of 3. Conversion to ionization rates (^) 
and energy deposition rates per unit of hydrogen density (Hx/n) 
is very sensitive to the adopted X-ray spectrum (Wardle et al., 
1999); adopting ^(s^') ~ 2 x 10"' (erg s^' cm^^-,^ j„fgr 
^ ~ (0.2 - 1) X 10"'^ S-' and Hx/n < (1 - 5) x 10"^° erg cm^ 
S-' (Maloney et al., 1996) for NGC 4418-Arp 220 in the far-IR 
nuclear photosphere traced by the OH. It is unclear if such low 
energy deposition rates per unit density can produce the high 
nuclear OH abundance we measure at least in NGC 4418, where 
very high densities are found. Nevertheless, cosmic rays are an 
alternative, and both X- and cosmic rays could efficiently gen- 
erate OH from a high reservoir of H2O that is formed by other 
means (i.e. by a hot core chemistry or by low-velocity outflows), 
through photodissociation of H2O by internally generated UV 
photons (Gredel et al., 1989; Wardle et al., 1999) and/or reac- 
tions of H2O with He"", H+, H+, and H+ to form OH directly or 
via the H20'^ and H^O^ intermediates; this high H2O reservoir is 
evident in both sources. Detailed chemical modeling is required 
to differentiate between these, and other alternatives; while H2O 
may be expected to be abundant in regions where HCN-NH3 or 
OH are enhanced, it remains to be determined whether the HCN- 
NH3 and the OH can be cospatial. 

4.2.3. Sources of luminosity 

The presence of high columns of HCN and H2O and the implied 
hot core chemistry in the nuclear regions neither supports nor 
rules out the presence of a deeply buried AGN. It could simply 
indicate that, regardless of the nature of the intrinsic luminos- 
ity source(s), frequency degradation in a very thick medium (as 
UV radiation is absorbed by dust and re-emitted at longer wave- 
lengths by samphng progressively colder media), makes far-IR 
photons more relevant for the chemistry of large fractions of the 
nuclear gas (through dust heating and mantle evaporation) than 
the quickly absorbed UV photons. 

Taking the Galactic Orion hot core as a template, its size 
(» 10" diameter, or 0.02 pc at 414 pc; e.g. Wilson et al., 2000) is 
~ 900 and > 4000 times smaller (and the square of these num- 
bers in surface) than those estimated for the Ccore and Cwest com- 
ponents of NGC 4418 and Arp 220, respectively. If the HCN ab- 
sorption is interpreted in terms of hot cores around young stars, 
then ~ 8 X lO"* and ~ 1.5 x 10^ cores are inferred, respectively. 
Likewise, if a typical hot core harbours a luminous star with 10^ 
Lq, ~ 7x10^ and ~ 7x10'' young stars are inferred, respectively. 
This interpretation, however, has the drawback that Galactic hot 
cores are short lived, with a time scale of order 10^ yr (WUner 
et al., 2001). Even if the high pressure in the nuclear regions 
(Spoon et al., 2005) allows the hot cores to survive longer, a 
synchronized birth of ~ 10'' hot cores in < 1 Myr over spa- 
tial scales of > 80 pc, involving a propagation velocity of stel- 



lar formation of > 100 kms , seems unlikely. The observed 
absorption is rather interpreted in terms of large hot-core like 
structures (disk/torus/envelope) not necessarily associated with 
the very first stages of star formation. 

The inferred hot core chemistry, complemented by the strong 
effects of X- and/or cosmic rays, is found to be physically asso- 
ciated with warm and luminous sources of far-IR emission. The 
surface brightnesses of the nuclear regions in both sources are 
very high, ~ 7 x lO'^ Lo/kpc^ over 20 pc for the Ccore compo- 
nent of NGC 4418, and (0.9-4) x 10'^ Lo/kpc^ over 150-80 pc 
for the western nucleus of Arp 220. In the latter source, these are 
lower limits if the dust is heated from inside (Downes & Eckart, 
2007). Sa8 indicates that the luminosity from the Arp 220 west- 
ern nucleus could be generated by 200-1000 super-star clusters. 
An important contribution from star formation to the luminos- 
ity of the western nucleus is indeed evident (e.g.. Smith et al., 
1998; Lonsdale et al., 2006), and the apparent lack of strong X- 
ray emission from X-ray binaries associated with star formation 
(Lonsdale et al., 2006) can be interpreted in terms of extreme 
absorption column densities (Sa08, this work). 

On the other hand, the broad 1.6 GHz OH megamaser emis- 
sion terminal velocity of 800 kms"' in Arp 220 (Baan et al., 
1989), comparable to those of OH in the far-IR absorption lines 
of AGN driven ULIRGs (Fischer et al., 2010; Sturm et al., 201 1), 
is much higher than the OH velocities observed in the far-IR 
spectrum of Arp 220. We interpret this as evidence for both 
far-IR extinction of the high velocity gas accelerated in the in- 
nermost nuclear region and deceleration of the outflowing gas 
with increasing involved columns. A very compact source for 
the outflowing gas would also support the presence of an ener- 
getically significant, strongly buried AGN that would also ex- 
plain the somewhat different distributions traced by the com- 
pact centimeter-wave sources and the submilhmeter continuum 
emission (Sa08). In NGC 4418, the extreme HCN and H2O ex- 
citation on a scale of 20 pc in diameter, the large fraction of 
luminosity coming from such a compact region, and the indi- 
cations of internal heating, favor the scenario of a far-IR thick 
disk/torus/envelope surrounding a nascent AGN, as has been 
proposed in previous works (e.g. Spoon et al., 2001), or a hyper- 
star cluster ~ 70 times more luminous than the super-star cluster 
in NGC 5253 (Turner et al., 2000). 

In NGC 44 1 8 , the outward acceleration due to radiation pres- 
sure can be described by Orad ~ kctT^/c ~ (2 - 5) x 10"'' cm 
s~^, where ^ ~ 2 - 5 cm^ g"' is the Rosseland mean opacity 
at 100 - 150 K (Thompson et al., 2005) and cr is the Stefan- 
Boltzmann constant. We find that a^d is similar to the gravi- 
tational acceleration, ag^v = GM/R^ ~ (2.5 - 5) x 10~* cm 
s~^ where M = (1.6 - 3.2) x 10^ Mo (for 7200 -1-2) and 
R = 10 pc have been estimated (Table 1). If the stellar mass is not 
much higher than the gas mass in the Ccore component, the above 
comparison suggests that radiation pressure plays an important 
role in supporting the nuclear structure against self-gravity in an 
Eddington-like limit (Scoville, 2004). In Arp 220 (Cwest) the im- 
portance of radiation pressure appears less important due to the 
higher M/R^ and lower* T^. Nevertheless, the disk in the western 
nucleus may be inclined such that our view of the warmer dust 
is obscured in the far-IR and we have a "rotationally- supported" 
edge-on view of the structure. 



Vertical radiation pressure support in optically thick nuclear disks 
imply kT^^ cc (e.g. Thompson et al., 2005), where T^fi is the effective 
dust temperature and 2^ is the gas column. Comparison of NGC 4418 
(Ccore) and Arp 220 (Cwest) indicates, however, lower T^s and higher Eg 
in Arp 220. 
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4.3. The extended component ofArp 220 

The extended and halo components (Cextended and Chaio) in 
Arp 220 reveal a chemistry very different from that found in the 
nuclear region. The H2O abundance drops to the relatively low 
value of ~ 5 X 10 typical of Galactic PDR regions, and the 
OH/H2O ratio reverses the nuclear value becoming > 1, as in 
Galactic translucent clouds (e.g. Neufeld et al., 2002). Our esti- 
mated OH abundances of several X 10^^ are in general agreement 
with those derived by Polehampton et al. (2005) in the intervin- 
ing translucent clouds toward Sgr B2. On the other hand, the out- 
flow signatures observed in the OH 163 and 79 jum doublet, in 
the [O i] 63 and 145 jjxa lines (Fig. 11), and in the H2O 3o3 -> 2i2 
line (Fig. 4p), appear to indicate that the outflow is spatiaUy ex- 
tended on the back side of the nuclei. 

The extended component of Arp 220 resembles the translu- 
cent Galactic medium, though with a higher ionization rate as 
derived from the absorption in lines of molecular ions, as we 
wiU report in a forthcoming paper. 

4.4. Inflowing gas In NGC 4418? 

The intringuing redshifted absorption in the ground-state lines of 
OH and O i, and the blueshifted emission in the [O i] 63 jum line, 
caused us to consider in §3.6 the possibility that an extended en- 
velope around the nucleus of NGC 4418 is inflowing onto the 
nuclear region, thus feeding the nucleus and providing an exter- 
nal ram pressure that may further confine the nuclear gas. Krug 
et al. (2010) have reported the detection of inflows in Seyferts 
based on absorption measurements of the Na I doublet at optical 
wavelengths. 

The morphology of the region as seen at near-lR wavelengths 
(Evans et al., 2003) and the absorption depths indeed suggest 
that the redshifted absorbing component is more than a single 
cloud in the galaxy disk with an anomalous velocity; the red- 
shifted absorption affects the whole column of gas external to 
the nuclear region (~ 10^^ cm"^), as otherwise the ground-state 
OH absorption at the systemic velocity would be deeper than 
observed. Thus the OH and O i absorptions suggest bulk inward 
gas motions. 

If spherical symmetry is assumed for this component, a mass 
inflow rate of 12 Mq yr"' is obtained at the surface of the nu- 
cleus, based on the model shown in Fig. 1 1 . Spherical symme- 
try is consistent with the OH line profiles, in particular with 
the emission in the OH 163 jum doublet. However, the spheri- 
cal approach overpredicts the [O i] 63 fim blueshifted emission 
(Fig. 11a). Furthermore, pure radiative transfer cannot rule out 
that a significant fraction of the [O i] 63 fim emission feature is 
formed in the nucleus, rather than in the envelope, if the high 
density component (several xlO^ cm"^, as derived from OH, 
H2O, and HCN) is sufficiently rich in atomic oxygen (a few 
xlO'^). Observational surveys of the [O i] 63 //m line in Galactic 
hot cores may help further evaluate the latter possibility, but pure 
hot core chemistry predicts O i to have an abundance of only 
< 10"^ (Charnley, 1997). Guided by these models, we do not 
expect much contribution from dense hot cores to the [O i] 63 
fim emission feature. However, the high-lying OH lines are trac- 
ing a region where photodissociation of H2O is probably taking 
place (§4.2.2), so OH photodissociation could produce signifi- 
cant columns of dense O i. StiU, the C-shock models by Wardle 
et al. (1999) with enhanced H2O photodissociation due to X- or 
cosmic rays produce little O i, with OH/Oi > 100, but models 
for XDRs and CRs by Meijerink & Spaans (2005) and Meijerink 



et al. (2011) predict large amounts of oxygen in atomic form, 
though they use densities significantly lower than 3x10'' cm"-'. 

In summary, our observations and analysis favor an inflow in 
NGC 4418, though the geometry and inflow rate of the inflow- 
ing gas are rather uncertain. If the gas were inflowing only from 
the front side of the nucleus with an opening angle of, say, 90°, 
the mass inflow rate would be of only « 2 yr"'. Kawara et 
al. (1990) first suggested a possible inflow in NGC 4418 based 
on the narrow CO / = 1 - line. The interaction with the 
nearby galaxy VV 655 could have trigered the inflow via a close 
(perigalactic), non-merging encounter ~ 10^ yr ago for an av- 
erage velocity difference of ~ 200 kms"' in the plane of the 
sky, and models (Noguchi, 1988; lono et al., 2004) predict that 
the inffow triggered by interacting galaxies may last that long. 
Alternatively, a non-axisymmetric gravitational potential as in 
barred galaxies can fuel the nuclear region with large amounts 
of gas (Friedli & Benz, 1993), which is supported by the pecu- 
har morphology of NGC 4418. 

4.5. Evolutionary effects 

In spite of the very different global parameters (mass, luminos- 
ity, size, and merger versus non-merger class) of Arp 220 and 
NGC 44 1 8 , their nuclear regions share common properties . They 
contain warm dust, very excited molecular absorption, and high 
surface column densities of H2O, OH, HCN, and NH3. There 
are, however, significant differences in the properties revealed 
by their far-infrared spectra that probably stem from and further 
probe their individual evolutionary histories. The smaller, less 
massive, and less luminous NGC 4418 shows warmer dust and 
higher molecular excitation than Arp 220. Perhaps the most sug- 
gestive spectroscopic differences between the two sources are 
their inferred kinematics and '*'0 to '**0 ratios: while NGC 4418 
appears to have a ratio similar (to within a factor of 2) to that 
seen in the local solar neighbourhood, the ratio in Arp 220 is 
significantly lower. 

The outffow signatures in Arp 220, which range from 800 
kms~^ in its OH megamaser radio emission to 100 kms~^ in 
submillimeter transitions of HCO^ and CO and up to 500 km s~' 
in some moderate excitation submillimeter and far-lR emission 
lines, suggest that injection of enriched gas from current and 
previous stellar processing has occurred. On the other hand, me- 
chanical feedback has not yet occurred in NGC 44 1 8, which may 
still be in a stage of inflow. Analysis of the H^^'O and the '^'OH 
fines indicate "'O-to-'^O ratios of 70 - 130 in Arp 220, remi- 
niscent of the low ratio of ~ 40 derived for Mrk 231 (Fischer 
et al., 2010). The strong absorption in the red A-component of 
the ground-state '^OH 120 jum doublet (Fig. 8c) indicates that 
the '^O enhancement is not restricted to the nuclear region, but 
also apphes to the Cextended and Chaio components. At milhme- 
ter wavelengths, Martm et al. (201 1) derive a ^^O-to-^^O ratio of 
> 80 in Arp 220, while Harrison et al. (1999) estimated a ratio of 
~ 150 in NGC 253. In NGC 4418, however, we find that flie '''O- 
to-i^O ratio is > 250, wifli a best fit value of ~ 500. While ^^O is 
a primary product of stellar nucleosynthesis, '^O is a secondary 
one, and this ratio decreases with increasing star forming gener- 
ations (Prantzos et al., 1996). Since galactic outflows are domi- 
nated by entrained ISM gas (Veilleux et al., 2005), the low ratios 
found in Arp 220 and Mrk 23 1 probably trace stellar processing 
throughout the merger lifetime. The apparent ISM inflow, lack 
of outflowing gas, and higher ^^O-to-^^O ratios in NGC 4418 
appear to indicate an earlier evolutionary phase in the Ufetime of 
the central AGN or nuclear starburst. 
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We tentatively suggest that the nuclei of Arp 220 and 
NGC 4418 are in an earlier evolutionary phase than the QSO 
driven ULIRG Mrk 231 (e.g. Veilleux et al., 2009). In this 
source, high velocity (~ 1000 kms"') feedback (Fischer et al., 
2010; Sturm et al., 2011; Rupke & Veilleux, 2011) in the form 
of a massive molecular/atonuc outflow has more strongly dis- 
rupted its nuclear intersteUar medium, as evidenced by the lower 
dust opacity of its warm component (Gonzalez- Alfonso et al., 
2010) compared with the nuclear regions of NGC 4418 and 
Arp 220. The lack of measurable HCN emission/absorption in its 
Herschel/SPIRE spectrum (van der Werf et al., 2010), suggests 
that exposure to UV and X-ray emission from its AGN appears 
to have at least partially shut down the hot core nitrogen chem- 
istry^. Further, analysis of Herschel/SPIRE observations of H2O 
emission in Mrk 23 1 (Gonzalez- Alfonso et al., 2010) suggest an 
H2O abundance of ~ 10"^, significantly lower than in NGC 4418 
and Arp 220, but observations of the same lines would help ver- 
ify this issue. Lastly, the stronger stellar processing suggested by 
the increased enhancement of "^OH in Mrk 231 may indicate an 
older age of its merger induced star-formation than in Arp 220. 
This evolutionary scenario is consistent with the fact that in the 
mid-IR Mrk 231 looks AGN-like, while Arp 220 and NGC 4418 
are still dominated by deep silicate absorption indicative of a 
buried, earher evolutionary stage (Spoon et al., 2007). 

5. Conclusions 

The sahent observational results presented here on NGC 4418 
and Arp 220 are: 

- Full range PACS spectroscopy of NGC 4418 and Arp 220 
reveals high excitation in H2O OH, HCN, and NH3. In 
NGC 4418, absorption fines were detected with Slower > 800 
K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), 
while in Arp 220 the excitation is somewhat lower with 
Slower > 600 K (H2O), 600 K (OH), 650 K (HCN), and 400 
K (NH3). 

- More specificafiy, NGC 4418 shows 38 H2O fines in ab- 
sorption, and with marginal detections of even higher-lying 
lines, while Arp 220 shows absorption in 28 H2O lines. All 
strong OH lines within the PACS domain are detected in both 
sources. In NGC 4418, at least 6 HCN rotational lines, from 
7 = 18 ^ 17 to y = 23 ^ 22 with Slower from 650 to 1075 
K are detected in absorption, while the 7 = 18 ^ 17 and 
possibly the / = 21 — > 20 lines are detected in Arp 220. 
Rotational NH3 lines arising from Slower up to ~ 600 K are 
detected in NGC 4418, while Arp 220 shows NH3 fines up 
to ~ 400 K. Molecular absorption equivalent widths from 
high-lying levels are higher in NGC 4418 than in Arp 220. 

- Up to 10 lines of Hj^'O and 3 '^OH A-doublets are detected 
in Arp 220, while 6 H^'^O fines are detected in NGC 4418, 
where only a marginal detection of ''^OH is found. The lines 
of the rare isotopologues are thus significantly stronger in 
Arp 220 tiian in NGC 4418. 

- Lines with moderate excitation in Arp 220 show indications 
of outflowing gas, specifically in the redshifted emission fea- 
tures observed in lines of H2O, OH, and O i. The absorption 
in more excited molecular lines peak at central velocities rel- 
ative to our adopted redshift of 0.0181, though a contribu- 
tion to the absorption by outflowing gas with velocities up 

' But not completely: during the final revision of this manuscript, 
Aalto et al. (2011) reported the detection of bright HC3N (10-9) in the 
nuclear region of Mrk 231. 



to ~ 150 kms"' is not ruled out. In NGC 4418, however, 
the molecular lines do not show any indication of outflowing 
gas, and the absorption in the low-lying OH doublets as well 
as in the [O i] 63 fiia line are redshifted. 

We have analyzed the observed absorption-donunated line 
profiles together with the far-IR continuum, using multicompo- 
nent radiative transfer modeling with a minimum number of pa- 
rameterized components. The conclusions stemnning from the 
analysis are: 

- Extremely high H2O column densities per unit of continuum 
optical depth at 50 jum (T50), > 3 x 10'** cm"^ and OH/H2O 
and HCN/H2O ratios of « 0.5 and 0.1 - 0.3, respectively, 
are derived in the nuclear region of NGC 4418, along with 
dust temperatures ~ 150 K in the very inner nuclear core. 
In the nuclear region of Arp 220, Tj is in the range 90 - 
130 K, implying H2O columns of (6 - 0.8) x 10'** cm^^ per 
T50 and OH/H2O and HCN/H2O ratios similar to those of 
NGC 4418. The HCN/NH3 ratios are 2 - 5 in both sources. 
The continuum emission in these nuclear regions is optically 
thick at far-IR wavelengths in both galaxies. 

- The nuclear molecular abundances relative to H are esti- 
mated to be high, with X112O ~ 10"^ in NGC 4418 and 
(0.2- 1.5) X 10-5 in Arp 220. 

- The absorption by the "^O isotopologues indicate a strong 
enhancement of '**0 in Arp 220, with a '^0-to-"^0 ratio of 
70- 130, while the ratio in NGC 4418 is significantiy fiigfier, 
> 250 with a best fit value of ~ 500. The apparent ^^O 
enhancement in Arp 220 probably traces steUar processing 
throughout the merger lifetime. 

- We also require an extended component to produce the ab- 
sorption in the low-lying fines of OH and H2O; the frac- 
tional H2O abundance in Arp 220 drops in this region to 
~ 5 X 10"^ and the OH/H2O is reversed relative to the nu- 
clear region to 2.5 - 10. In the envelope of NGC 4418, we 
find;^'OH - a few x 10"^. 

- The extreme H2O, HCN, and NH3 columns and abundances 
in the nuclear region of both galaxies, together with the high 
dust temperatures, are interpreted in terms of a hot core 
chemistry, where grain mantles are evaporated. The high OH 
columns and abundances require in addition photoprocesses 
due to X-rays and/or enhancement by high cosmic ray fluxes. 
H2O and OH may also be boosted in low-velocity shocks 
with enhanced H2O photodissociation. 

- Very thick (A^h ^ 10^^ cm"^) nuclear media are responsible 
for the frequency degradation of the intrinsic luminosity in 
both sources, as high energy photons are absorbed and re- 
emitted at far-IR wavelengths sampling progressively colder 
media. Such conditions, and not necessarily a starburst- 
dominated luminosity, may also be responsible for the hot 
core chemistry in the nuclear regions. The high surface 
brightnesses, ^ 6x 10'^ and (0.9 - 4) x lO'^ Lo/kpc^ on 
spatial scales of ~ 20 and 150 - 80 pc in NGC 4418 and 
Arp 220, respectively, indicate either strongly buried AGNs 
or hyper star clusters. In NGC 4418, radiation pressure is 
probably responsible for the support of the nuclear structure 
against self-gravity. 

- Both the outflowing gas and the low '^^0-to-'**0 ratio in Arp 
220 suggest a stage more evolved than in NGC 4418, but less 
evolved flian flie QSO driven ULIRG Mrk 23 1 . 

- In NGC 4418, the redshifted absorption observed in the low- 
lying OH and [O i] lines is best explained in terms of an 
envelope inflowing onto the nuclear region, but the geometry 
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of this inflow is rather uncertain; we estimate the mass inflow 

rate to be < 12M0yr"'. 
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Table 3. H2O line equivalent widths, fluxes, and velocity shifts. 



Line 




A 
(jim) 


(kms"-') 


NGC 4418 

Flux" 
(10-2' Wcm-2) 


Khift 
(km s-') 


(fans'-') 


Arp220 

Flux" 
(10-2' Wcm-2) 


(km s"') 


H2O 533 


-^422 


53.138 


24.9(6.5) 


-20.0(5.2) 


35(29) 


1.9(2.6) 


-4.5(6.3) 




H2O 827 


^7,6= 


55.131 


6.2(2.3) 


-4.8(1.8) 


-3(16) 








H2O43, 


— >322 


56.325 


35.6(1.8) 


-26.6(1.4) 


-8(33) 


30.0(1.2) 


-65.3(2.6) 


-29(49) 


H2O 422 


->3l3 


57.636 


37.1(2.1) 


-28.0(1.5) 


-15(35) 


22.8(1.3) 


-49.3(2.8) 


14(43) 


H2O 432 


^321 


58.699 


50.7(0.9) 


-37.8(0.7) 


16(34) 


34.7(1.0) 


-74.9(2.1) 


-87(47) 


H2O726 


^615^ 


59.987 


5.6(1.3) 


-4.1(1.0) 


-6(22) 








H2O431 


^4o4 


61.809 


11.1(2.6) 


-7.3(1.7) 


37(24) 


2.3(1.4) 


-4.9(3.0) 




H2O 818 


— >7o7 


63.324 


13.3(1.0) 


-9.1(0.7) 


33(18) 


2.4(0.6) 


-5.1(1.3) 




H2O 7,6 


-^625 


66.093 


11.8(0.9) 


-7.5(0.6) 


-6(24) 


5.7(1.0) 


-12.3(2.1) 




H2O 330 


^221 


66.438 


64.7(1.3) 


-41.8(0.9) 


-18(34) 


40.8(0.8) 


-87.8(1.8) 


-55(38) 


H20 33, 


— >220 


67.089 


43.6(2.2) 


-28.1(1.4) 


12(31) 


41.6(0.9) 


-90.4(2.1) 


-30(49) 


H20 330 


— > 3o3 


67.269 


13.2(1.3) 


-8.5(0.8) 


-9(29) 


7.7(2.3) 


-16.7(4.9) 


-24(51) 


H2O 524 


^4i3 


71.067 


23.9(1.7) 


-13.6(1.0) 


-45(24) 


16.7(1.0) 


-33.2(2.0) 


35(44) 


H2O717 


^6o6 


71.540 


17.6(1.8) 


-9.7(1.0) 


27(34) 


9.5(0.8) 


-19.0(1.6) 




H2O 7o7 


^6i6 


71.947 


25.5(2.3) 


-13.8(1.2) 


-11(32) 


10.0(0.8) 


-20.1(1.7) 


35(38) 


H20725 


^634'= 


74.945 


5.0(1.5) 


-2.6(0.8) 


-15(23) 








H2O 321 


^2i2 


75.381 


73.9(1.1) 


-38.6(0.6) 


-23(32) 


79.7(0.8) 


-153.4(1.6) 


4(41) 


H2O 423 


^3,2 


78.742 


40.3(1.9) 


-19.1(0.9) 


-36(27) 


34.4(0.9) 


-62.9(1.6) 


-42(43) 


H2O 615 


— >524 


78.928 


19.4(1.8) 


-9.3(0.9) 


-4(27) 


4.9(0.8) 


-8.8(1.4) 


-70(31) 


H2O 616 


^5o5 


82.031 


31.4(1.4) 


-14.1(0.6) 


2(25) 


20.3(0.6) 


-35.6(1.1) 


-21(37) 


H2O 606 


^5i5 


83.284 


25.1(1.7) 


-11.3(0.7) 


-16(27) 


16.1(1.2) 


-27.5(2.0) 


-4(41) 


822 


^2ii 


89.988 


60.4(1.7) 


-23.3(0.7) 


5(27) 


39.8(1.1) 


-60.6(1.6) 


-32(35) 


H2O 643 


^634' 


92.811 


8.0(1.6) 


-3.0(0.6) 


-35(19) 








H2O 625 


^6i6 


94.644 


10.1(1.7) 


-3.6(0.6) 


-37(17) 








H2O441 


^432 


94.705 


14.3(1.8) 


-5.1(0.6) 


-33(18) 








H2O 5l5 


^4o4 


95.627 


42.9(2.4) 


-15.2(0.8) 


-11(30) 


22.8(2.1) 


-31.6(2.9) 


9(41) 


H2O 440 


^431 


95.885 


12.5(2.0) 


-4.3(0.7) 


-2(20) 


3.0(1.7) 


-4.1(2.3) 




H20 6,s 


^605" 


103.940 


10.0(1.4) 


-2.9(0.4) 


68(29) 


7.6(1.4) 


-9.5(1.7) 


-34(42) 


H2O 22, 


^lio 


108.073 


67.4(2.6) 


-17.5(0.7) 


-23(46) 


72.2(2.1) 


-82.9(2.4) 


40(47) 


H2O 524 


^5,5 


111.628 


24.8(2.2) 


-6.1(0.5) 


42(55) 


4.2(0.7) 


-4.5(0.8) 




H20414 


^3o3 


113.537 


36.7(2.1) 


-8.5(0.5) 


14(39) 


26.0(2.9) 


-27.0(3.1) 


-32(44) 


H2O404 


^3i3 


125.354 


25.4(1.8) 


-4.6(0.3) 


-19(33) 


27.7(1.2) 


-23.7(1.0) 


-6(54) 


H2O 331 


^322" 


126.714 


12.0(1.8) 


-2.1(0.3) 


55(35) 


8.8(5.3) 


-7.4(4.4) 


14(119) 


H2O 423 


^4,4 


132.408 


28.1(2.5) 


-4.6(0.4) 


5(68) 


12.6(0.8) 


-9.7(0.6) 


-29(46) 


H20514 


— >5o5 


134.935 


26.1(1.9) 


-4.1(0.3) 


0(49) 


8.6(0.8) 


-6.3(0.6) 


18(40) 


H20 330 


-^321 


136.496 


36.9(1.9) 


-5.5(0.3) 


0(46) 


19.8(0.8) 


-14.3(0.6) 


-37(48) 


H2O 3 13 


— >2o2 


138.528 


30.7(1.4) 


-4.5(0.2) 


30(39) 


29.0(1.2) 


-20.2(0.8) 


-40(45) 


H20413 


^322 


144.518 


11.9(1.7) 


-1.5(0.2) 


-32(28) 


11.6(0.9) 


-7.3(0.6) 


-24(40) 


H20431 


^422 


146.923 


12.2(2.1) 


-1.5(0.3) 


-17(33) 


3.3(0.6) 


-2.0(0.3) 


-23(27) 


H2O 322 


^3l3 


156.194 


9.2(2.3) 


-1.0(0.2) 


30(26) 


14.5(0.9) 


-7.5(0.5) 


-29(38) 


H2O 3o3 


^2i2 


174.626 


13.2(2.5) 


-1.0(0.2) 


-7(15) 


20.3(3.7) 


-7.3(1.3) 


-33(26) 


H2O 2i2 


^loi'' 


179.527 


81.9(4.9) 


-5.5(0.3) 


118(35) 


167.3(5.1) 


-51.5(1.6) 


34(40) 



" Numbers in parenthesis indicate 1 cr uncertainties from Gaussian fits to the lines. 
Velocity shifts are relative to z = 0.00705 for NGC 4418 and z = 0.0181 for Arp 220. 
Marginal detection. 
Probably contaminated. 
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Table 4. [O i] and OH line equivalent widths, fluxes, and velocity shifts. 















NGC 4418 






Arp 220 




Line 








A 




Flux" 


,,a,b 
V 

shift 


eq 


Flux" 


wa.b 
shift 










(fim) 


(Km s ) 


/in— 21 wT 2\ 
(lU W cm ) 


(Km s ) 


(Km s ) 


/in— 21 m ^.^— 2\ 
(lU W cm ) 


(Km S ) 


[O I] 1 - 


2 






63.184 


25.2(1.0) 


-17.2(0.7) 


115(16) 


45.4(2.0) 


-96.3(4.3) 


20(30) 


[0 i]0- 


1 






145.525 


-7.3(1.6) 


0.9(0.2) 


35(22) 


-25.4(1.1) 


15.6(0.7) 


75(77) 


OHn3/2 


- 1^3/2 


5 " 
2 , 


3 + 
2 


119.233 


58.4(2.0) 


-12.0(0.4) 


125(36) 


212.5(2.2) 


-204.5(2.1) 


41(52) 


OH Hj/z 


- n3/2 


5 + 
2 


3" 
1 + 


119.441 


54.9(1.9) 


-11.2(0.4) 


114(32) 


174.9(2.0) 


-168.1(1.9) 


39(44) 


OHni/2 


-Hj/z 


1 - 

2 . 


3 + 
2 


79.118 


17.8(1.4) 


-8.7(0.7) 


97(24) 


63.5(1.3) 


-115.0(2.3) 


90(35) 


OH ni/2 


- n3/2 


1 + 

?4. 


3- 

2 


79.181 


12.5(1.3) 


-6.1(0.6) 


72(19) 


63.5(1.3) 


-115.0(2.3) 


90(35) 


OHHi/z 


- n3/2 


3 + 

2 


3- 

?-t 


53.261 


38.8(4.5) 


-31.3(3.6) 


76(35) 


95.7(3.8) 


-230.5(9.1) 


66(44) 


OHni/2 


- n3/2 


3- 
2 


3 + 

2 


53.351 


59.2(5.6) 


-47.7(4.5) 


56(52) 


94.4(3.7) 


-227.3(8.8) 


31(42) 


OHni/2 


- n3/2 


3 

2 , 


96.340 


9.1(3.3) 


-3.1(1.1) 


38(34) 








OHn3/2 


- n3/2 


7 + 
2 




84.420 


77.8(1.9) 


-34.2(0.8) 


10(30) 


88.6(1.4) 


-150.8(2.3) 


2(38) 


OHn3/2 


- n3/2 


7- 
2 


5 + 
2 , 


84.597 


89.6(2.0) 


-39.2(0.9) 


1(34) 


95.4(1.4) 


-161.7(2.3) 


-16(39) 


OHn3/2 


- n3/2 


1 + 


7 + 
2 


65.132 


74.1(1.7) 


-48.7(1.1) 


16(39) 


49.8(1.2) 


-107.7(2.6) 


-6(50) 


OH n3/2 


- n3/2 


9 ^ 
2 


7 

5? 
2 


65.279 


62.7(1.8) 


-40.8(1.1) 


-13(39) 


45.8(1.2) 


-98.8(2.5) 


-57(49) 


OHni/2 


-ni/2 


7 

!i" 


71.197 


72.1(2.8) 


-40.5(1.6) 


-16(45) 


42.4(1.0) 


-84.3(2.0) 


-17(54) 


OH n3/2 


- n3/2 


2 


9- 


52.934 


24.1(5.0) 


-18.8(3.9) 


-57(22) 


3.0(2.2) 


-7.4(5.5) 


-67(40) 


OH n3/2 


- n3/2 


11 - 

2 


7 0^ 


53.057 


23.5(5.4) 


-18.5(4.3) 


91(25) 


3.6(3.8) 


-8.8(9.3) 


4(33) 


OHni/2 


-ni/2 


9 

!" 


2 


55.920 


54.1(4.4) 


-41.0(3.3) 


40(91) 


12.5(2.5) 


-26.9(5.4) 


-52(58) 


OHHi/z 


-ni/2 


2 


1 - 


163.124 


-10.0(2.2) 


1.0(0.2) 


-78(24) 


-71.9(2.7) 


31.9(1.2) 


150(77) 


OH ni/2 


-ni/2 


3- 
2 


2 


163.397 


-9.6(2.2) 


0.9(0.2) 


-16(23) 


-41.8(5.0) 


18.5(2.2) 


161(57) 



" Numbers in parenthesis indicate 1 <t uncertainties from Gaussian fits to the lines. 
" Velocity shifts are relative to z = 0.00705 for NGC 4418 and z = 0.0181 for Arp 220. 
The two A-components are (nearly) blended into a single spectral feature. 



Table 5. HCN line equivalent widths, fluxes, and velocity shifts. 



Line 






A 
(jim) 


W" 
(km s-') 


NGC 4418 
Flux" 
(10-2' Wcm-2) 


ya,b 
shift 

(km s-') 


W" 
(km s'') 


Arp 220 
Flux" 
(10-21 Wcm-2) 


ya,b 
' shift 

(km s^') 


HCN 25 


— » 


24 


135.626 


5.8(1.8) 


-0.9(0.3) 


-20(43) 








HCN 23 




22 


147.365 


9.4(1.9) 


-1.1(0.2) 


69(23) 








HCN 22 


— » 


21 


154.037 


12.3(1.6) 


-1.3(0.2) 


19(35) 


4.3(2.0) 


-2.3(1.1) 


-35(58) 


HCN 21 




20' 


161.346 


14.7(2.7) 


-1.5(0.3) 


7(21) 


9.6(1.5) 


-4.5(0.7) 


-53(45) 


HCN 20 




19 


169.386 


12.7(2.7) 


-1.0(0.2) 


-6(20) 








HCN 19 




18 


178.274 


11.0(3.8) 


-0.8(0.3) 


21(27) 


3.8(2.0) 


-1.2(0.6) 


24(72) 


HCN 18 


— » 


17 


188.151 


21.0(5.4) 


-0.9(0.2) 


-16(20) 


15.2(3.2) 


-3.1(0.6) 


-100(37) 



" Numbers in parenthesis indicate 1 cr uncertainties from Gaussian fits to the lines. 
" Velocity shifts are relative to z = 0.00705 for NGC 4418 and z = 0.0181 for Arp 220. 
" Doubtful identification in Arp 220. 
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Table 6. NH3 line equivalent widths and fluxes for relatively isolated spectral features. 











NGC 4418 




Arp220 


Line 








(km s"') 


Flux" 


(km s-') 


Flux" 








(A'tn) 


nO"-' W cyvr^~\ 


(10-2' Wcm-2) 


NH3 (8, 6)i- 


^ (7, 6)a 




62.888 


3.9(1.7) 


-2.7(1.1) 






NH3 (7, 6)a 


^ (6,6),v 




72.439 


6.7(0.6) 


-3.6(0.3) 


4.6(0.5) 


-9.2(1.0) 


NH3 (6, 5)s 


^ (5, 5)a 




83.432 


10.1(1.2) 


-4.5(0.5) 


4.5(0.8) 


-7.6(1.4) 


NH3 (6, 4)s 


^ (5, 4)a 




83.590 


13.6(1.4) 


-6.1(0.6) 






NH3 (6, 3)s 


^ (5, 3)a 




83.709 


10.1(1.0) 


-4.5(0.5) 


3.0(0.7) 


-5.1(1.1) 




^ (5, 2)a 




83.794 


4.5(1.0) 


-2.0(0.4) 


b 


b 


NH3 (6, 1 - 


0)i ^ (5, 1 - 0)a'^ 




83.844-83.861 


11.0(1.1) 


-4.9(0.5) 


7.5(1.2) 


-12.9(2.1) 


NH3 (4, 3)i 


^ (3,3)a 




124.648 


17.9(1.7) 


-3.3(0.3) 


25.1(1.1) 


-21.6(0.9) 


NH3 (4, 2 - 


1-0)5-^(3,2-1- 


0)a^ 


124.796- 124.913 


38.3(1.6) 


-7.0(0.3) 


21.5(1.2) 


-18.5(1.1) 


NH3 (4, 3 - 


2 - l)a-> (3,3-2- 


1)5= 


127.108 - 127.181 


38.9(1.3) 


-6.9(0.2) 


28.9(0.9) 


-24.1(0.7) 


NH3 (3, 2)i' 


^ (2, 2)a'i 




165.597 


30.9(2.3) 


-2.8(0.2) 


27.6(2.0) 


-11.8(0.9) 


NH3 (3, \)s 


^ (2, 1)0" 




165.729 


13.6(2.4) 


-1.2(0.2) 


b 


b 


NH3 (3, 2 - 


1-0)0^(2,2-1- 


0)5^ 


169.968 - 169.996 


26.9(3.4) 


-2.2(0.3) 


33.7(1.5) 


-13.2(0.6) 



" Numbers in parenthesis indicate 1 cr uncertainties from Gaussian fits to the Unes. 
Uncertain due to baseline subtraction or blending. 
Joint emission of several A^-lines 
Probably contaminated by H3O*. 



Table 7. Hj^O line equivalent widths and fluxes. 



Line X 
(jum) 



H'*0 432^32i 59.350 

331 ^220 67.883 

H^80 32i^2i2 75.867 

h™ 423 ^3 12 79.531 

H^*'0 322^2i, 90.936 

H^^'O 22,^1 10 109.347 

Hf0 4i4^3o3 114.296 

H'«0 3i3^2o2 139.586 

h|80 3o3^2i2 174.374 

Hf0 22i^2i2 183.530 



NGC 4418 

Flux" 

(kms-') (10-"Wcm-2) 



1.5(0.9) -1.0(0.6) 

2.0(0.9) -1.2(0.6) 

7.0(1.6) -3.5(0.8) 

6.0(1.1) -2.9(0.5) 

3.5(1.5) -1.3(0.6) 

5.7(1.6) -1.5(0.4) 

b b 



Arp220 



"eq 

(km s-') 


(10-21 W cm-2) 


7.9(1.1) 


-17.2(2.5) 


8.3(0.9) 


-18.0(1.9) 


17.9(0.8) 


-34.7(1.5) 


5.4(0.8) 


-9.6(1.4) 


4.2(0.7) 


-6.4(1.1) 


13.2(1.8) 


-14.8(2.0) 


7.0(1.5) 


-7.2(1.5) 


10.0(0.8) 


-6.8(0.5) 


7.9(1.5) 


-2.9(0.6) 


14.8(2.5) 


-3.8(0.7) 



" Numbers in parenthesis indicate 1 cr uncertainties from Gaussian fits to the lines. 
* Uncertain due to baseline subtraction. 



Table 8. '*OH line equivalent widths and fluxes. 



Line 








X 
(jum) 


NGC 44 18 

Flux" 

(km s"') (10-21 W cm-2) 


W" 
(km s-i) 


Arp220 

Flux" 
(10-21 W cm-2) 


''OHn3/2 


- 113/2 


9 + 
2 


7- 
2 


65.690 


1.7(0.9) 


-1.1(0.6) 


4.9(0.9) 


-10.6(1.9) 


'*'OHn3/2 


- 113/2 


7 + 
2 




84.947 


2.1(0.9) 


-0.9(0.4) 


7.0(0.8) 


-11.8(1.3) 


'*OH n3/2 


- n,/. 


7- _ 


h 


85.123 


0.7(0.7) 


-0.3(0.3) 


12.9(0.8) 


-21.5(1.4) 




- Ui/2 


^ + 




120.171 


0.4(0.8) 


-0.1(0.2) 


38.7(2.2) 


-36.7(2.0) 



Numbers in parenthesis indicate 1 cr uncertainties from Gaussian fits to the Unes. 



